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ABSTRACT 

The extraordinary 1998 August 27 giant flare places strong constraints 
on the physical properties of its source, SGR 1900+14. We make detailed 
comparisons of the published data with the magnetar model, which identifies the 
Soft Gamma Repeaters as neutron stars endowed with ~ 10^^ G magnetic fields. 
The giant flare evolved through three stages, whose radiative mechanisms we 
address in turn. The extreme peak luminosity L > lO^Ledd, hard spectrum, and 
rapid variability of the initial ~ 0.5 s spike emission all point to an expanding 
pair fireball with very low baryon contamination. We argue that this energy 
must have been deposited directly through shearing and reconnection of a 
magnetar-strength external magnetic field. Low-order torsional oscillations of 
the star fail to transmit energy rapidly enough to the exterior, if the surface 
field is much weaker. A triggering mechanism is proposed, whereby a helical 
distortion of the core magnetic field induces large-scale fracturing in the crust 
and a twisting deformation of the crust and exterior magnetic field. After the 



^Canadian Institute for Theoretical Astrophysics, 60 St. George St., Toronto, ON M5S 3H8 
^University of Texas, Department of Astronomy & McDonald Observatory, Austin, TX 78712, USA 
■^Institute for Theoretical Physics, University of California, Santa Barbara, CA 93106, USA 



- 2 - 



initial spike (whose ~ 0.4 s duration can be related to the Alfven-crossing 
time of the core), very hot (T ^ 1 MeV) plasma rich in electron-positron pairs 
remains confined close to the star on closed magnetic field lines. The envelope 
of the August 27 flare can be accurately flt, after ~ 40 s, by the contracting 
surface of such a "trapped fireball." The form of this fit gives evidence that the 
temperature of the trapped pair plasma decreases outward from its center. We 
quantify the effects of direct neutrino-pair emission on the X-ray light curve, 
thereby deducing a lower bound /irnin ~ 1 x 10^^ G-cm^ to the magnetic moment 
of the confining field - comparable to the strongest fields measured in radio 
pulsars. The radiative flux during the intermediate ~ 40 s of the burst appears 
to exceed the trapped flreball flt. The lack of strong rotational modulation and 
intermediate hardness of this smooth tail arc consistent with the emission from 
an extended pair corona, in which 0-modc photons are heated by Compton 
scattering. This feature could represent residual seismic activity within the 
star, and accounts for ~ 10 percent of the total flare flucnce. We consider in 
detail the critical luminosity, below which a stable balance can be maintained 
between heating and radiative cooling in a conflned, magnetized pair plasma; 
but above which the conflned plasma runs away to a trapped flreball in local 
thermodynamic equilibrium. The emergence of large- amplitude pulsations at 
~ 40 s probably represents a transition to a pair-depleted photosphere whose 
main source of opacity is electrons (and ions) ablated from the heated neutron 
star surface. The best flt temperature of the black body component of the 
spectrum equilibrates at a value which agrees well with the regulating effect 
of photon splitting. The remarkable four-peaked substructure within each 
5.16-s pulse, and the corresponding coUimation of the X-ray flux, has a simple 
explanation based on the strong inequality between the scattering cross sections 
of the two photon polarization modes. The width of each X-ray 'jet' is directly 
related to the amount of matter advected outward by the high-cross section 
ordinary mode. 

Subject Headings: gamma rays: bursts - stars: neutron - X-rays: stars 



1. Introduction 

On 1998 August 27, a giant flare of hard X-rays, lasting more than 5 minutes, was 
detected by Konus-Wind, Ulysses and BeppoSAX (Hurley et al. 1999a; Feroci et al. 1999; 
Mazets et al. 1999). It was only the second giant outburst detected from the Soft Gamma 
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Repeaters (SGRs), following nearly 20 years after the famous event on 1979 March 5 (Chne 
et al. 1982). The four or five known SGRs are identified by their frequent emission of 
much less energetic and shorter-duration bursts {E < 10^^ ergs, At <, 1 s), whose peak 
luminosities reach nevertheless ten thousand times the Eddington luminosity. Statistically, 
the short bursts are similar to earthquakes and solar flares (Cheng et al. 1996; Gogus et 
al. 1999; Gogus ct al. 2000). 

SGR 1900+14, the source of the August 27 outburst, was first detected in 1979 
(Mazets ct al. 1979a), and became active again in 1992 (Kouveliotou et al. 1993) and in 
1998 May (Kouvehotou et al. 1998a, Hurley et al. 1999b). ASCA and RXTE observations 
after the 1998 May activity episode revealed a periodicity of 5.16 s in the quiescent X-ray 
emission with period derivative ~ 1 x 10~^°ss""^ (Hurley et al. 1999c, Kouvehotou et 
al. 1999; Woods et al. 1999). A rapid spindown P ~ 10"^° s s"^ and 7.47 spin period 
was also measured for SGR 1806-20 (Kouveliotou et al. 1998b) These apparent rotation 
periods are similar to that of SGR 0526-66, which displayed an 8 s periodicity during the 
giant flare of 1979 March 5 (Mazets et al. 1979b, Barat ct al. 1979). At least two of the 
SGRs are associated with supernova remnants (SNRs) and/or radio plerions (Hurley 2000). 
However, the recently- verified position of SGR 1900+14 (Vasisht et al. 1994; Hurley et 
al. 1999d; Murakami et al. 1999; Frail, Kulkarni & Bloom 1999) lies just outside the edge 
of G43.8+0.6, a < 10^-year-old galactic SNR. It is plausible that SGR 1900+14 was born 
in this SNR, since the other known giant flare source, SGR 0526-66, has a position inside 
but near the the edge of a <, 10^-year-old SNR in the Large Magellanic Cloud (Cline et 
al. 1982). 

The SGRs have been identifled with magnetars (Duncan & Thompson 1992, "DT92"): 
neutron stars in which magnetism rather than rotation or accretion is the dominant source 
of energy for radiative and particle emissions. The hypothesized ~ 10^^ — 10^^ G magnetic 
flelds of these stars are consistent with many observed SGR properties, including the very 
high energies and hyper- Eddington luminosities of the outbursts, the long spin periods, 
rapid spindown rates, strong X-ray emissions in the quiescent state, and the young ages 
of active stars (DT92; Paczyhski 1992; Thompson & Duncan 1995, "TD95"; Thompson & 
Duncan 1996, "TD96"; Thompson et al. 2000). 

A recent joint paper (Feroci et al. 2000, hereafter Paper I) analyzed the Ulysses and 
BeppoSAX observations of the August 27 giant fiare, and drew some conclusions about 
the underlying physics of the event. The intense, hard spike of gamma-rays during the 
first 0.4 second was identified with a freely expanding pair fireball, which requires a very 
clean energy source. We found that the subsequent rotation-cycle averaged flux diminished 
monotonically, or nearly so. This strongly suggests that most of the burst energy was 



-4- 



released at early times, near the beginning of the flare, and that the energy which did not 
escape promptly during the initial hard spike (t < OA s) was retained by the star as a 
residue, to leak away on longer timescales. 

In principle, such an energy residue could take several forms. Some energy could be 
in the form of an optically-thick pair-photon plasma contained in the magnetosphere by 
closed magnetic field lines (TD95), which we will call a "trapped fireball." Additional 
energy could be stored within the star itself in the form of a internal oscillations, including 
a torsional Alfven mode in the liquid core; or a standing shear wave in the neutron star 
crust (McDermott, Van Horn, & Hansen 1988; Duncan 1998). Thus, observations of the 
flare and its aftermath give clues about the way in which the energy residue is partitioned 
among various components. 

In this paper, we begin by discussing the general issue of SGR energy budgets in §0, 
estimating the minimum magnetic field needed to power both the bursting and quiescent 
emission of the SGRs. Then, in we consider the physical states of magnetars on the 
threshold of a giant flare, with a particular focus on how and where the energy that powers 
the flare is stored, and how its catastrophic release may be triggered. Our favored model 
involves a wound-up magnetic field that is pinned by the neutron star crust. 

With this background information developed, we then begin to consider successive 
stages of the 1998 August 27 flare. In Paper I we showed that the envelope of the August 
27 light curve, smoothed over the 5.16-s rotation period during the interval from 40 to 400 
s after the onset of the flare, is well-fit by a mathematical curve describing the cooling of a 
zone of very hot plasma composed of 7-rays and electron-positron pairs. Such a "trapped 
fireball" cools through the formation of a sharp temperature gradient inside its outer 
boundary, which propagates inward in a "cooling wave" (TD95). In particular, this model 
matches the rapid, late decline in the X-ray flux at ~ 400 s after the flare onset. The 
observed terminal drop in the X-ray flux is does not appear to be consistent with models 
in which the emissions are powered by cooling surface hotspots or by incremental internal 
mode damping. In ^ of this paper, we relate the measured slope of the lightcurve to the 
internal properties of the fireball and study the effects of neutrino cooling on the trapped 
fireball. 

The extreme properties of the giant flare lead to stringent bounds on its magnetic field 



(DT92; Paczyhski 1992; TD95). In §5.1| we use the condition that the trapped fireball 
energy must remain confined by the magnetic field to derive lower bounds on the star's 
magnetic dipole moment. We generalize this argument by considering deviations from 
centered-dipole geometry, which imply a lower magnetic moment for a fixed fireball energy. 
We show that the lower bound on the magnetic moment cannot be significantly weakened 
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without raising the temperature of the confined plasma sufficiently that almost all its energy 
is lost to neutrino pair emission, rather than X-rays. 

Then, in ^5 .21 , we show how the extremely high peak luminosity of the hard spike (more 
than a million times the Eddington luminosity of a neutron star) points to the presence 
of magnetic fields stronger than ~ 10^'* G. Because the magnetic field lines are tied to 
material of a very high density, a sudden release of elastic stresses in the deep crust occurs 
only gradually compared with the time ~ -Rns/c ~ 3 x 10""^ s for an Alfven-like excitation 
to cross the magnetosphere. As a result, the high measured luminosity is consistent with 
external shearing and reconnection of the external field, but only if the fiux density exceeds 
~ 10^^ G. For essentially the same reason, an internal shear wave cannot transmit energy 
to the magnetosphere at ~ 10^ Lcdd in weaker magnetic fields. We constrast our results 
with the previous suggestion that an internal p-mode or /-mode excitation of a neutron 
star could induce strong shock heating of the surface layers and electromagnetic damping 
in the magnetosphere (Ramaty et al. 1980; Lindblom & Detweiler 1983), even if the surface 
magnetic field were only ~ 10^^ G. We explain how two basic assumptions of this model - 
that the surface layers would be shock heated, and that external Alfven modes could be 
excited with 6B/B ~ 1 by an internal oscillation - are very likely incorrect. Indeed, Blaes 
et al. (1989) have shown that internal shear modes have a much stronger coupling to the 
magnetosphere than do p- and /-modes. 

Some aspects of giant fiare emission spectra and spectral evolution are elucidated in 
§1^. In particular, we consider the hard spectrum of the intial spike in § |6.1| , and argue that 
the hardest 7-ray photons must be generated at more than ~ 10^" cm from the star, due 
to the same limitations of pair opacity that arise in models of cosmological gamma-ray 
bursts. The equilibrium temperature of the pulsating soft tail fitted in Paper I is shown 
in § |6.2| to agree with the calculation in TD95 of the freeze-out of photon splitting in a 
strongly-magnetized scattering atmosphere. 

The first ~ 40 seconds of the August 27 fiare showed only mild and irregular pulsations 
at the 5.16 spin period, with a higher fiux and harder spectrum than the remainder of 
the burst. This excess emission over the predictions of the trapped fireball curve amounts 
to 20% of the burst fiuence (Paper I). In §^ we interpret this 40-s "smooth tail" as the 
signature of a transient and extended pair corona, powered by residual creep of the neutron 
star crust following the principal disruption. In particular, we consider the thermodynamic 
stability of a strongly magnetic pair atmosphere, in which X-ray photons are subject 
simultaneously to Compton scattering and splitting. We show that this corona is thermally 
unstable if the pairs are heated at a rate exceeding ~ 10^^ ergs s~^ within a volume 
(10 km)3. 
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Around ~ 40 s following the flare trigger, large-amplitude modulations appeared in the 
light curve with a remarkable four-peaked structure. These peaks were maintained with a 
5.16-s period (equal to the spin period measured in quiescence) during most of the rest of 
the event. In we examine in more detail the proposal of TD95 that the cooling X-ray 
flux from a trapped flreball in a super-QED magnetic fleld will be coUimated along (partly) 
open magnetic fleld lines. In this picture, the complicated pulse shape provides a template 
of of higher multipoles in the neutron star's surface magnetic fleld. 

Finally, in §^ we summarize all of the evidence, based on the study of giant flares and 
their aftermaths, for magnetic flelds ^ 10 -Bqed = 4.4 x 10^'* G in SGRs, and the main new 
implications of this work for the physics of magnetars. 

In what follows, we use the convenient reference fleld 

Sqed = ^ = 4.4 X 10^=^ G, (1) 

at which the non-relativistic Landau energy heB/rrieC formally becomes equal to the 
electron rest energy rrieC^. 



2. Strong Magnetic Fields in the SGR Sources 

The electromagnetic output of an SGR source consists of both hard X-ray outbursts 
and the persistent X-ray emission. The core of the star may release an even larger energy in 
the form of thermal neutrinos, over its active liftime (TD96). We now estimate the minimal 
magnetic fleld needed to power the observed activity. The close affinity between the SGR 
and AXP sources (some of which have variable persistent X-ray output even in the absence 
of bursts) indicates that an SGR can dissipate magnetic energy in two distinct modes: first, 
in an episodic manner through brittle fracturing of the crust; and, second, through a more 
plastic deformation of the star (TD96). 



2.1. Energetic Requirements 

The association of several magnetar candidates (three AXPs and at least two SGRs) 
with young supernova remnants, tsNR ^ 10"^ yr (Hurley 2000), combined with the short 
spindown times of these stars, suggests an active lifetime of order of magnitude ~ 10^ yr. 
Two giant fiares have been detected from four verified, active SGR sources within the last 
30 years. The implied fiaring rate is ~ 10~^ per year per source, suggesting a total number 
of giant fiares A^^giant ~ 10^ from each source during its active lifetime ~ 10'^ yr. Since the 
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measured X-ray and gamma ray output of the 1979 March 5 and 1998 August 27 flares 
exceed ~ 5 x 10^*^ ergs and ~ 1 x 10^^ ergs respectively, the total flare energy released per 
source is 

^-^3X10«(^^) (^) ergs (2) 

The persistent X-ray luminosity of SGR 0526-66 is ~ 10^^ ergs s~^ (Kulkarni et al. 
2000) and is approximately lO^^ ergs s"^ for SGR 1806-20 and SGR 1900+14 during periods 
of quiescence (Murakami et al. 1994; Woods et al. 2000). The net energy released in 
persistent X-ray emission can therefore be estimated as 

Ex - Lx tactivc = 3 X 10^^ f ^) ergs (3) 

yiO"*^ ergs s yiO^yry 

over the active lifetime. Thus, the observed electromagnetic output of an SGR is fairly 
evenly divided between flares and its highly non-thermal persistent emission. The 
observation of a transient enhancement in the persistent output of SGR 1900+14 following 
the August 27 flare (Woods et al. 2000) establishes a close physical connection between the 
two phenomena. This enhanced emission has been ascribed to persistent currents excited by 
non-potential deformations of the external magnetic fleld, caused by the shifting positions 
of magnetic footpoints which are anchored to an evolving crust (Thompson et al. 2000). 

It is possible that signiflcant power is emitted at wavelengths intermediate between 
X-rays and optical in sources such as SGR 0526-66 with soft power-law spectral components 
(Kulkarni et al. 2000). Such intermediate wavelengths are not observable in known galactic 
magnetar candidates because of the opacity of the intervening interstellar medium. In the 
case of the low-extinction LMC source SGR 0526-66, spectral measurements at <, 0.5 keV 
are not yet available. For sources with hard power-law spectra (photon index ~ —2 or less 
steep) extending to low energies, the correction factor to eqn. (^ would be of order unity 
only. 

It is instructive to compare -Eflare with the energy detected in ordinary, short-duration 
SGR bursts. These events are distributed over a wide range of energies ~ 10^'^ — 10"^^ ergs, 
with a power-law dM jdE ~ E'^/"^^ so that the integrated energy is dominated by large 
events (Cheng et al. 1996; Gogus et al. 2000). Nonetheless, many fewer than ~ 10^ events 
with energies ~ 10^^ ergs have been detected during the same time that two giant flares 
have been observed. This suggests that the time-averaged output in short bursts is lower 
than either in the giant flares or persistent emission. 
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2.2. Magnetic Field Strength 

For the reasons just described, the persistent and bursting output of an SGR source 
appear to draw from the same energy reservoir, which at minimum is E^[^ ~ 10^^ ergs per 
source. The present rotational energy is much smaller, = 5 x 10^^(P/6 s)~^ ergs. It has 
been suggested that the SGR/AXP sources were born spinning much more rapidly, P ^3 
msec, and that a fraction of this initial rotational energy was converted to a strong magnetic 
field through an a — f2 dynamo in the convective proto- neutron star (DT92). This magnetic 
field would then provide a reservoir of energy for later bursting activity. An alternative 
form of potential energy for bursting activity involves an elastic deformation of the crust, 
resulting from the changing figure of the rotating star as its rate of spin decreases (e.g Baym 
et al. 1969). However, the maximum elastic energy which can be stored by the crustal 
lattice is only '^M^mstV^^l = 1 x 10^^ ^y^/^QS g-i)2 (M^mst/W-^ Mq) (^,,/0.003)2 
ergs (eq. []TT|), too small to power a single giant flare except under the most optimistic 



assumptions about the breaking strain ipcr- (Here is the shear wave velocity, which is 
nearly independent of density in the deep crust.) Even in sources which have not emitted 
giant flares, the elastic energy would need to be replenished on a timescale short compared 
with the spindown age; but the present rotational energy is insufficient for that purpose. 

Accretion power is not a promising energy source for the hard-spectrum, hyper- 
Eddington SGR outbursts, nor for persistent X-ray emission from SGRs. The observation 
of persistent emission from SGR 1900+14 within 1000 s of the 27 August 1997 giant 
flare (Woods et al. 2000) seems inconsistent with accretion-powered scenarios, because 
even ~ 10~^ of the radiative momentum of the flare would excavate the accretion flow 
and suppress accretion for a much longer time period (Thompson et al. 2000). Accretion 
models also have difficulty accomodating the power-law distribution of energies of ordinary 
SGR bursts, and the lognormal distribution of waiting times between bursts (Cheng et 
al. 1996; Gogus et al. 1999, 2000). Recent HST bounds on the optical emission from 
SGR 0525-66 set an upper limit of L^pt/Lx ~ 10~^ for reprocessed optical emission from 
an accretion disk (Kaplan et al. 2001), an order of magnitude less than is expected from 
an accretion disk (Perna, Hernquist, & Narayan 2000). Optical and infrared observations 
also exclude the simplest disk models of the Anomalous X-ray Pulsars [AXPs] (Hulleman 
et al. 2000a; Hulleman, van Kerkwijk & Kulkarni 2000), which resemble SGRs in inactive 
phases. In particular, 4U 0142+61 has a possible optical counterpart with a much redder 
spectrum than is typical of LMXBs (Hulleman et al. 2000b). Some other arguments against 
accretion-powered models of SGRs are given in §7.3 of TD95. 



Magnetic energy, on the other hand, can be a good "clean" power source for both 
hyper-Eddington bursts and the hard fireball component of the giant flares (DT92; 
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Paczynski 1992; TD95). A first estimate of tlie minimum magnetic field is obtained by 
averaging over the volume of the star, 

/ 6Ejam \ ^ o inl4 -1/2 / -^min \ ^ f -^NS \ '^^^ , , 

S^in = = 8 X 10^^ TT^ir-— 7771— G. (4) 



.eBUfisJ ^ V 10^^ ergs/ VlOkm 

Here denotes the efficiency of conversion of magnetic energy to X-rays. 

If the protons in the stellar interior form a Type II superconductor, then the magnetic 
field is confined to fluxoids with flux density comparable to the lower critical field 
Bel ~ 2 X 10^^ G (e.g. Easson & Pethick 1979). The mean field in the core is related to 
Bel by the volume- filling fraction / of the fluxoids, (B) = Bdf, and the magnetic energy 
per unit volume is B"^^/ /8n = {B)Bci/8n. When the mean field (B) < B^i (that is, when 
the efficiency of conversion of core magnetic energy to X-rays is less than ~ 10 percent) 
expression (^) decreases by a factor ~ e^^"^ {B^^^/ Bd). Otherwise, the minimum field is 
not modified to first order by proton superconductivity. 

If eB ~ 1, expression (^) is similar to the dipolar magnetic fields that are inferred 
from the rapid spindown of the SGRs (Kouveliotou et al. 1998b, 1999) and AXPs (TD96). 
However, the total magnetic energy of the star is almost certainly larger than its minimum 
value, because the currents that support the field flow through very highly conducting 
material and cannot entirely dissipate. Very strong core magnetic fields also seem to be 
required by models of core ambipolar diffusion which include the accelerating effects of core 
heating by the decaying field (TD96; Heyl & Kulkarni 1998). Only fields stronger than 
~ 10^ -Bqed ~ 5 X 10^^ G are transported over a distance ~ 10 km in times less than ~ 10"^ 
yr, typical of SGR active lifetimes. 

The magnetic flux density inside an SGR/AXP source probably lies below the upper 
critical field where proton superconductivity is fully quenched, Bc2 ~ 10^^ G. Nonetheless, 
the extended magnetic sheaths of the fluxoids are very densely packed in fields ;^ i?ci- As a 
result, field transport driven by collective effects such as vortex-line sweeping (Ruderman 
1991) is probably suppressed in slowly-rotating magnetars. 

Several hues of evidence suggest that SGR bursts and flares are powered by internal 
magnetic stresses acting on the deep crust and core. The hard initial 7-ray spikes observed 
in the two giant flares had durations ~ 0.2 — 0.5 s. This is comparable to the time for 
a >j 10^^ G magnetic field to rearrange material in the dense stellar core (TD95), but it 
is much longer than the timescales for global rearrangements of either the magnetosphere 
(~ -Rns/c ~ 0.03 ms) or the stellar crust (-Rns/^ ~ 10 ms; see also eqs. and pl[ 



below). Furthermore, the cumulative fluence of short SGR bursts is observed to rise linearly 
with time during periods of activity (Palmer 1999). This "relaxation system" behavior 
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indicates that a reservoir of energy is steadily loaded and undergoes stochastic, catastrophic 
relaxations. For example, such a reservoir could encompass a fraction of the star's surface, 
steadily dragged from below by the tension of field lines experiencing core ambipolar 
diffusion (TD96). This model is roughly analagous to what happens in earthquakes, where 
relaxation-system behavior driven by crustal plate motion is also found (see Palmer 1999; 
Gogus et al. 1999, 2000, and references therein). 

The magnetic field geometry of an SGR/AXP source almost certainly involves higher 
multipoles (for which the August 27 lightcurve provides clear evidence: Paper I; below); 
and quite plausibly a strong toroidal field component in the deep crust and core, as would 
result from post-collapse a-Q dynamo action (DT92; Thompson & Duncan 1993, hereafter 
TD93). This suggests that the magnetic energy of a magnetar is dominated by the internal 
field, with a probable value 



where Jb is the fraction of the star's volume filled by the strong field. 

Comparisons of eqs. (^, (H) and (^ indicate that only a fraction es ~ 0.03//^ of the 
star's magnetic energy must be dissipated in the form of X-rays and gamma-rays to power 
all observed SGR activity. Indeed, most of the core field energy is probably lost to thermal 
neutrinos or ultimately remains trapped in the core when the epoch of SGR activity ends. 



The crust of a magnetar is subjected to strong, evolving magnetic stresses. The star 
must initially relax to an equilibrium configuration before the crust forms (but after it 
becomes compositionally stratified: Lattimer & Mazurek 1981; Reisenegger & Goldreich 
1992). Observations of magnetic white dwarfs prove that static, magnetized equilibrium 
states exist even in the absence of any rigidity due to solidification. After the neutron 
star crust forms, the crust, core and field evolve through a sequence of equilibrium states 
in which magnetic stresses are balanced by both hydrostatic forces and elastic stresses in 
the crust. The star evolves via the very slow transport of the field through the core (by 
ambipolar diffusion) and the crust (by Hall drift) (Goldreich & Reisenegger 1992; TD96). 
Because the crust has a finite shear strength, these equilibria are punctuated by starquakes 
whenever the crust is strained past its breaking point. 




ergs. 



(5) 



3. 



Giant Flare Mechanism 



We now review the forms in which potential energy can be stored inside a magnetar, 
and how some of this stored energy can be transferred to the external magnetic field. 



to heat, and to internal torsional and standing shear waves. This suggests a candidate 
mechanism for triggering giant flares like the March 5 and August 27 events. 



3.1. Stored Elastic and Magnetic Energy 

A Soft Gamma Repeater stores a certain amount of potential energy which can 
drive rapid rearrangements of the external magnetic field, thereby triggering bright X-ray 
outbursts. This potential energy can be divided into an elastic component in the crust, and 
magnetic components in the crust, core and magnetosphere. The portion of the magnetic 
energy which is available to do mechanical work can easily be estimated, in a situation 
where the accumulation and release of stress is caused by pinning of the internal magnetic 
field by the crust. 

The elastic energy corresponding to a static shear strain ifj can be well approximated 

^-elas = 1 fltfj' r^dr = i (^^] (^') Mcoulomb V^,. (6) 

2 J j_Rns-a_Ro z y7 — (J.z J 

In this expression, the shear modulus isQ 

/i = 1.1 X 10^° p^J ergs cm-^ (7) 

and the shear wave speed, 

Vf. = {p/pf^ = 1.0 X 10=^pr4°-'km s-\ (8) 

depends weakly on the density p = pu x 10^^ g cm~^ at densities above neutron drip 
(Baym & Pines 1971; Strohmayer et al. 1991, Fig. 3). The integral @ extends from the 
base of the Coulomb lattice at a depth ARq (density po and shear wave speed V^o)- It has 
been evaluated using the equation of hydrostatic equilibrium dP/dr = —pg, and fitting a 
power-law relation between the pressure and density profiles, P(r) oc [p(r)]'^. In particular, 
the equation of state of Negele & Vautherin (1973) gives 7 ~ 1.47 and 

P(p) = 3.75 X 10^2 plf ergs cm'^ (9) 



"'This generalizes the expression given in Duncan (1998), to ahow for a variable position (density) at the 
base of the rigid Coulomb lattice. 

^The neutron star crust is here approximated as a body-centered cubic (bcc) Coulomb lattice (Ogata & 
Ichimaru 1990; Strohmayer et al. 1991), with ionic mass and charge determined by a power-law fit to the 
equation of state of Negele & Vautherin (1973). 
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above a density 3.6 x 10 gm cm . The mass within the Coulomb lattice 

Mo„„,„„. = r- p(r)i.rHr = ^.^3sEM (10) 

excludes a transition layer at its base where the nuclei have strongly aspherical shapes, 
forming rod-like and plane-like structures, with bulk elastic properties resembling those of a 
liquid crystal (Pethick & Potekhin 1998). Normalizing the base density po to a fraction 0.6 
of the nuclear saturation density psat = 2.7 x 10^^ g cm~^, expressions (D and (|10D become 

\ilO-^)y \0.Q Ps.t) VlOkmy ^1-4^0^ ^' ^ ^ 

and 



Mc„„,„„.^ 0.025A-^^J Afe. (12) 

The scaling with stellar radius i?NS and mass M^s assumes a fixed equation of state, and 
arises because the hydrostatic pressure -P(po) ~ [fl'(-RNs)^Couiomb/4vri?^g] at the base of the 
crust is fixed by the nuclear interactions (Lorenz, Ravenhall, & Pethick 1993). The factor 

= 1 — 2GMns/-Rnsc^ — 0.6 — 0.7 is a relativistic correction to the surface gravity, 
g{R^s) = GMnsA/Z^^s- 

Comparing expression with the energy of the 1979 March 5th event, 
-^Marchs ~ 5 X 10"^"^ ergs, it is clear that even for the optimistic case ipcr ~ 10~^, 
giant flares cannot be powered by pure elastic energy stored within the crust of a neutron 
star. Furthermore, it is often assumed that tpcr is one or more orders of magnitude less than 
the value ~ 10^^ appropriate for a perfect bcc lattice, since lattice imperfections inevitably 
weaken the solid (e.g. Ruderman 1991). But the critical strain ipcr for a large-scale yield 
in the gravitationally-stratified crust of a neutron star is actually highly uncertain. Any 
motion of the crust over a scale of kilometers is severely constrained, because the hydrostatic 
pressure in the deep crust exceeds the shear modulus by a factor ~ 10^. Degeneracy 
pressure and buoyancy forces provide strong resistance to bulk compressions and to vertical 
displacements within the crust. Observations of giant flares probe the behavior of Coulomb 
solids in a regime of high pressure and large applied stress that has no direct experimental 
analogs. 

It should be emphasized at this point the magnetic field which stresses the crust can 
store much more potential energy than the Coulomb lattice itself. Within the crust, the 
force balance 

^ ~ ^p (13) 
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implies that 



when B < B^. Here 



{6Bf/87r _ f B^ 



-"crust 



(14) 



5,, = (47r/i)V2 = 4xlOiV?4' G (15) 

is a characteristic field above which the crust tends to respond plastically to applied 
magnetic stresses. For B > B^, the equilibrium condition (|I3]) requires that the strain 
angle of the Coulomb solid exceed 6B/B, even while the field remains tied to the highly 
conducting solid (TD95). In weaker fields B < 5^, the field stores more potential energy 
than the lattice. This field energy might be tapped in short SGR bursts, which would 
require that the crustal material undergoes abrupt plastic deformation. However, a very 
weak field B ^ B^ does not have the strength to move the crust, even as the ratio in 
eq. (p^ formally becomes enormous. For example, when 

B < = i^H' B, = lx 10^^ p[t G (16) 

magnetic stresses within the crust are too weak to induce any yields or fractures. 



3.1.1. A Globally Twisted Magnetic Field 

The giant fiares involve a large disturbance which probably is driven by a rearrangement 
of magnetic field in the deep crust and core (TD95). It is important therefore to consider 
the case in which the crust is stressed from below by the evolving magnetic field in the liquid 
interior of the neutron star. A related configuration has been considered by Ruderman 
(1991), in a situation where mean core field is much weaker than 10^^ G, and its transport is 
driven by spindown. As we now show, a factor similar to (|T^) relates the available magnetic 
energy stored in the core to the elastic energy built up in the crust. The fiux density B 
which enters this relation turns out to be the poloidal field which threads the core-crust 
boundary. 

At least two lines of argument suggest that the interior magnetic field of the SGR and 
AXP sources is strongly wound, with a large toroidal component: first, their association 
with neutron stars that were formed with rapid rotation and strong differential rotation 
(DT92); and, second, their non-thermal and transient persistent X-ray emission (Thompson, 
Lyutikov, & Kulkarni 2001). We therefore focus on a cylindrically symmetric star with a 
uniform poloidal field B^ < B^ in its interior (Fig. |I]). In the core, the evolving field is 
assumed to be twisted about the magnetic dipole axis (the axis of symmetry), generating 
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a field component B^. In static equilibrium, before the critical point is reached, the 
stress BzB^/Ati applied to the lower base of the crust is balanced by an elastic stress 
= fi^{d(f)/dw) = jjiip within the crust, so that 

- 'wa^A. (17) 



Att vj dw 

Here, Ai?^ is the stress-averaged depth of the crust, 

0.47 



M 



^ /i(po)^o \dpl pg ^ V0-6Psat; VlOkmy \IAMq^ 

(18) 

(Duncan 1998). Expression ( [iTl) remains a good approximation as long as the spherical 
curvature of the neutron star crust can be neglected; that is, at cylindrical radius w ^ |-Rns- 
Integrating out from the axis of symmetry (tu = 0), and assuming the power-law dependence 

137^0 oc '?/' oc ro", (19) 

equation (p!7|) implies B^B^/An = (1 + a){AR^/zu)pip. The parameter a = 2 corresponds 
to a twist angle that is independent of w; whereas a = corresponds to a current that is 
localized on the symmetry axis. 

Notice that the equilibrium strain ipi^) increases away from the axis of symmetry for 
a > (eq. ^). This means that the crust can be expected to break first at a radius zofj-ac 
where 

V'(^frac) ~ V'cr, (20) 

the critical strain angle. If the strain increases monotonically with time, this means that 
fracture first occurs at Wfrac ~ (| — 1)-Rns- 

The energy stored in the twisted magnetic field is 

5E^ag = (2i?NS ■ TTca') 1^ = vr(a + (^) {AR^fR^s, (21) 



\B, 



per logarithm of cylindrical radius IaJ. This can be expressed as 

. 4 X 10« (p^) f^^)" r^)%rg= (22) 



for a = 2. It is also worth comparing eq. ( p2D with the crustal strain energy. The available 
magnetic energy in the core is the larger by a factor 



(BI/Stt) i?NS ■ vrti72 ^5^x 2 (1 + a)^Rns^R^ 



[lljj'^fi)AR^-7lw'^ \Bj 



(23) 
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This expression exceeds (|T^) by a geometrical factor, which works out to ~ 2 (tJ7/|i?Ns)^^ 
in the case of a uniformly twisted field (a = 2). Thus, in the giant flare the main role of the 
crust is to serve as a gate to facilitate the storage and episodic release of interior magnetic 
energy, not as a reservoir of energy itself. 

In this toy model, the available magnetic energy (^2]) suffices to power a single giant 
flare, if ipcr ~ 10^^. This yield strain is larger than usually adopted for neutron star crusts 
(e.g. Ruderman 1991), but it may not be unreasonable given the large hydrostatic pressure 
and near incompressibility of the crustal material. The required value of ijjcr is reduced if 
the poloidal current flowing through the star varies in sign, with multiple interior zones in 
which the toroidal fleld has a magnitude comparable to its limiting value Then the 

winding energy ( [2^ ) is multiplied by the number of zones A^, and ipcr decreases by a factor 
N^^^"^ . However, the hydromagnetic stability of such a conflguration is not certain. 

Although the energy that can be released in one giant flare is limited by the shear 
strength of the crust, this simple model signiflcantly underestimates the total energy that 
can be released over successive flares. In addition, in this simplifled geometry the twist of 
the core fleld must be assumed to increase with time to reach a point of instability. Both 
apparent difficulties can be addressed by noting that the stratification of the core allows the 
field to retain a significantly larger twist than can be dissipated in a single fiare. 

Consider an electrically conducting medium that is stratified along planes, which 
run perpendicular to the direction of gravity. This stratification is assumed to be stable 
to convection, in the sense that internal g-modes have real frequencies. The gradient of 
electron fraction in the core of a neutron star provides just such a stable stratification 
(Reisenegger & Goldreich 1992). A cylindrical bundle of magnetic flux threads the medium, 
tilted at an angle a with respect to gravity (Fig. |2^). We assume that the magnetic fleld 
is weak, in the sense that the fluid motions driven by the J x B force lie almost parallel 
to the equipotential surfaces. If a = 0, then the flux bundle can unwind entirely through 
such constrained motions, as long as its ends are not pinned. By contrast, if a = 7r/2 and 
the flux bundle runs perpendicular to gravity, then any motion of the fleld and entrained 
conducting fluid which reduces the twist must involve the motion of fluid elements across 
the equipotential surfaces. This means (Fig. pE| ) that a localized twist will not be able to 
spread out along the flux bundle - even if the ends of the bundle are not pinned by external 
forces. 

For an arbitrary value of a, it is clear that any initial twist will be able to relax only 
partially. In particular, if the fleld is strongly twisted, ^ B^/ tana, then the flux bundle 
will relax to a magnetostatic equilibrium without B^/Bz being signiflcantly reduced in 
magnitude. To see this, consider a cylindrical bundle of flux which is unpinned at either 
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end. Coordinate z runs along the axis of the flux bundle, and gravity lies in the plane 
= 0. The twisted bundle is allowed to deform from this initial condition, subject to the 
constraint that the component of the velocity parallel to gravity vanishes. 



(We neglect deviations from cylindrical symmetry, so that in this coordinate system the 
radial components Bji and Vr both vanish. In effect, the flux bundle and the stratifled 
medium are assumed to be conflned with an inflnitely rigid cylindrical shell.) The induction 
equation then becomes 



One sees that if ^ B^j tana, then a rotation of surfaces of constant z cannot reduce 
B^ by more than a factor of two, because the second term in the induction equation is not 
symmetric about the axis of the cylinder. 

Slower transport processes such as ambipolar diffusion can change a over timcscales 
much longer than the hydromagnctic time -Rns/^Az- In this manner, clastic stresses can 
build up if the flux bundle is tied to a rigid medium (e.g. the neutron star crust). 



The two giant flares were initiated by hard 7-ray spikes whose ~ 0.2 — 0.5 s width is 
comparable to the Alfven crossing time of the star, if the internal poloidal magnetic fleld 
is ~ 10^^ G (TD95). By contrast, a relaxation of the external magnetic fleld has a much 
shorter characteristic timescale i?Ns/c ~ 3 x 10^^ s. This gives evidence that the giant 
flares are driven by the relaxation of internal magnetic stresses. 

Any motion of the crust of a neutron star is strongly constrained by degeneracy 
and buoyancy forces, which resist bulk compressions and vertical displacements. A 10^^ 
G magnetic fleld contributes only ~ 10~^ of the hydrostatic pressure at the base of the 
crust. At the same time, it is capable of deforming the crust along equipotential surfaces 
of the star. Over large scales (kilometers or more), the lowest energy deformations of the 
crust involve a displacement fleld ^ which has a vanishing radial component and satisfles 
V ■ (pO — 0) while maintaining a non-zero curl V x (p^) ^ 0. 

Such a large-scale twisting motion of the crust, if it occurs suddenly, probably involves 
the formation of one or more propagating fractures. This type of deformation has been 
associated with the giant flares (Thompson et al. 2000; TD95). Since the magnetic potential 



Vz cos a — V^ sin a cos = 0. 



(24) 




(25) 



3.2. Release of Stored Magnetic Energy 
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energy of twisting exceeds the crustal elastic energy (eq. [Q), the crust acts as a gate for 
the release of magnetic energy. The external magnetic field lines are anchored to the crust, 
so when the fracture occurs, the rapid turning motion does work on the external magnetic 
field. At the same time, it brings into contact regions of the crust where the magnetic field 
has differing strength and orientation, thereby creating strong field gradients and localized 
current sheets. Within the star, the bulk rotational motion also excites a helical Alfven 
wave in the liquid core (where the magnetic field provides the dominant restoring force to a 
torsional motion). Torsional shear waves are excited in the rigid crust (McDermott et al. 
1988; Duncan 1998), but are likely to carry less energy. We consider each of these physical 
elements in turn. 



3.2.1. External Shear and Reconnection 

Above the surface of the star, a disturbance of the magnetosphere propagates at very 
nearly the speed of light, which is some 300 times the shear wave speed in the deep 
crust. Thus, the external magnetic field can respond adiabatically to a smooth deformation 
of the crust. However, regions of the magnetosphere whose footpoints are strongly sheared 
can become subject to a purely magnetohydrodynamic instability (Lynden-Bell & Boily 
1994; Mijic & Linker 1994). If localized near a fault, such an instability probably leads to 
reconnection and induces magnetohydrodynamic waves outside the star. Transverse Alfven 
waves will have a characteristic frequency u; ~ c/x at a distance x from the fault. The 
corresponding wavenumbers are k_i_ ^ k\\ ~ a;~^ in directions parallel and perpendicular 
to the background magnetic field. These waves can damp rapidly by cascading to a high 
wavenumber through non-linear interactions (Thompson & Blaes 1998). When the rate 
of transfer of wave energy (the 'cascade luminosity') exceeds a critical value (~ lO^Lgdd 
within a volume of ~ (10 km)^), the dissipated wave energy is locked onto the magnetic 
field lines in a thermal pair-photon plasma (TD95; Section |^ below). If excited directly in 
the magnetosphere, these waves are not easily reabsorbed by the crust, since crustal shear 
waves of the same frequency have a much larger wavenumber, by the ratio c/V^ ~ 300. 

Thus, even though the giant flares are probably driven by internal magnetic stresses, 
the creation of a hot fireball outside the star involves the rapid relaxation of external 
magnetic stresses. This suggests that the energy released promptly in X-rays and 7-rays is 
roughly proportional to the external magnetic energy, after taking into account geometrical 
factors. Since a pure dipole field with polar strength Bi^ x 10^^ G carries an energy 
~ 10'^'' 5^5 (-Rns/10 km)^ ergs, the energy dissipated is small compared to the total external 
field energy if 5 ~ 10^^ G. 
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Suppose, for example, that a cap of the neutron star crust of radius w were to rotate 
through an angle A0, while the surrounding crust remains almost stationary. Even if this 
rotation were to be aligned with the axis of a purely dipolar field, one would expect strong 
shear to build up in the external field near the boundary of the cap. The amount of energy 
dissipated externally can be estimated, by supposing that this shear relaxes in N steps 
after a rotation through an angle A(f)/N. The magnetic energy released in each step is 
then ~ 7riu^(i?Ns/3) (5-BNs)^/8vr, where 5i?Ns/-BNS ~ A0/A^ and we have assumed a dipole 
geometry with a radial scale length i?Ns/3. The net energy which is available for powering 
X-rays and 7-rays is times this, 

^^-^^ Wl^j i?NS^NS-4xlO [j^) [j^) ergs. 

(26) 

Thus the fraction of the star's surface area that experiences slippage in a giant flare could be 
as small as (7rro^)/(47ri?Ns) ~ 0.02(iu/3 km)^ under favorable assumptions. If the position 
of the active region is determined by the large-scale winding of the internal magnetic field, 
then subsequent flares may occur at the same spot, and some of the surface may never 
flare. Alternatively, it is possible that independent regions encompassing no more than 
-^flare ~ 10^^ of the surfacc area flare only once over the active lifetime of an SGR. The peak 
luminosity of the initial spike is easily reproduced in this model, as long as exceeds 
the observed value. This leads to a lower bound on the surface field, as discussed in (§|5l^). 

In a more realistic case, high-order multipoles of the field are present. This implies the 
existence of some exterior magnetic field lines which remain close to the star in a complex 
geometry. Stretching of field lines whose footpoints lie on opposite sides of a fault provides 
another mechanism for fast reconnection and flaring. 

3.2.2. Internal Torsional and Standing Shear Waves 

A twisting deformation of the neutron star crust, driven by the Maxwell stress 
BzB^/Att, will be accompanied by a partial unwinding of the core magnetic field. We 
consider, as before, a uniform poloidal flux density Bz threading the core. Then the core 
supports a torsional Alfven mode, which propagates along the poloidal field with a period 

" = ^-HloI^J i l0^^gcm-3 j liol^j ^ (27) 

near the symmetry axis. Here, Vaz = Bz/y/4:7rp is the poloidal Alfven speed. This period 
receives a small correction from the shear strains that build up in the crust as the field that 
threads the crust/core boundary is twisted from below. 
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It should be noted that the mild density gradient in the neutron star core introduces a 
gradient in the wave period, with respect to the cylindrical radius w. This gradient has the 
effect of washing out any quasi-periodicity in the wave motion - and thence in the dynamic 
stress applied to the crust - after a dozen or so wave periods. 

Whether such a torsional Alfven mode is strongly excited during a giant flare depends 
on the rate at which the crust deforms or breaks. Let us suppose that the core field untwists 
through an angle A0 < A^q, where the net twist across the core is 

A^„(.) ^^El^^J^), ,28) 



This formula can be expressed as 



O.O2M0 / V108 cm s-V V10"V^3km/ VlO km. 

(29) 

in radians, making use of eq. ([T7|). If the crust moves slowly, over a timescale At much 
larger than the Alfven crossing time 2R^a,/VAz) then the amplitude of the resulting torsional 
wave is small. We estimate A0 

wave 

(A0/At) X (2i?Ns/VA2), or equivalently 

A0„ave 2_RnS /op|\ 

On the other hand, if the motion of the crust is slowed only slightly by friction 
(compared with a free torsional oscillation of the core) then a significant fraction of the 
released core magnetic energy goes into magnetic torsion modes. The fact that the observed 
durations of hard spikes in giant flares are comparable to a magnetar's torsional oscillation 
period, ~ 2i?Ng/VA2 = |-PAifven for A0 ~ A0o, is consistent with the hypothesis that 
the crust broke, with only modest frictional resistance after the event onset. Significant 
torsional excitation of the neutron star core is thus plausible: ii^wave ~ Ex ~ 10^'' — 10"^^ 
ergs. 

A similar conclusion holds for the proportion of the crustal strain energy (|ll|) which is 
converted to a torsional shear wave. The fundamental torsional mode (with n = radial 
nodes) has an amplitude which varies weakly with depth in the crust (McDermott et al. 
1988). Its harmonic is restricted to £ > 2 from angular momentum conservation, and its 
period is 
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(This period is shifted downward shghtly in the presence of a poloidal magnetic field, but 
the shift is not significant unless the component of the field along the direction of shear is 
comparable to 5^ ~ 6 x 10^^ G; Duncan 1998.) Notice that the period (^Tj) is significantly 
shorter than the torsional Alfven period (P?]). As a result, a significant fraction of the 
crustal strain energy would be converted directly to a torsional shear wave during a giant 
flare, only if this strain energy were released on a relatively short timescale compared with 
the duration of the hard spike, 0.03 s or less. Otherwise, the dynamic strain excited in the 
crust is dominated by the coupling to a torsional Alfven mode in the core. In either case, 
the amplitudes of the shear wave and Alfven wave can be expected to equilibrate rapidly 
(over a few wave periods) according to eq. ([13|), with the effect that the crustal shear wave 
has the lower energy by the factor (p^). 



3.2.3. External Torsional Deformation of the Magnetic Field 

Such an internal torsional mode, with frequency u = 27r/Pmode (eq. |^ or |^), also 
drives a torsional excitation of the magnetic field outside the star. In the case of the crustal 
shear mode, we focus on the fundamental mode, whose amplitude varies weakly with depth 
below the surface of the neutron star (McDermott et al. 1988). The surface amplitude of 
the external magnetic twist is related to the surface displacement ^ of the crust through 

^ = ^. (32) 

-DNS C 

We now show that, in equilibrium, only a tiny fraction of the energy of the internal mode 
is transferred to the external mode. This energy transfer could take two forms: first, an 
Alfven-hke mode propagating freely on very extended field lines; and, second, a quasi-static 
torsional deformation of the field lines closer to the star, which oscillates in sign at the same 
frequency as the internal mode. Field lines which can support a propagating Alfven mode 
have a minimum length cPmode, and extend out to a large radius -Rmax ^ |cPniode (in a 
dipole geometry). This works out to -Rmax ~ 10^° {Bz/lO^^ G)~^ cm, using the period (p^). 
These field lines are concentrated within a surface polar angle ^^(i?Ns) ^ 3i?Ns/cPmode- 

The internal mode will, in general, not be aligned with the external dipole moment, 
and so the surface shear will increase the external field energy by a fractional amount 
~ (^/i?Ns)^- This works out to 6E^^g{R > Pns) ~ i^^ns^ns (C/^ns)^ (here Pns is the 
polar dipole field). Comparing with the elastic energy i?eias = |'^^(^^)^mode stored in the 



crustal lattice, one finds 



^-£'mag(-R > -Rns) _ -^NS-^NS _5^]^q-8 ( -^NS \^ (^Rns\ ^ ( ^modc ^ ( -RnS 



^eias QM^oAe^^Rls G J \ c J \0.02MqJ VlOkm. 

(33) 

This ratio is small in the case of a crustal shear mode (Mmode — ^Couiomb = 0.02 M©) and 
even smaller in the case of a torsional Alfven mode (Mmode — ^ns)- 

The equilibrium energy in the freely propagating Alfven mode can easily be shown 
to be smaller than eq. (|33|) by a factor ~ Rm/Rma.x = 3-RNs/c-Pmode- Along the 
excited bundle of field lines, the product of the wave energy density {6BY/8n and the 
cross-sectional area A of the field lines is a constant, and proportional to the square 
of the current I{9) flowing at polar angles smaller than 9 ~ (-RNsZ-Rmax)^''^- One has 
I{9) = (c/47r) 7r6'^i?Ns(V ^ B)/j, where the current density is related directly to the surface 
shear through (V x B)^;; ~ 27r (^/i?Ns) -Bns/-Rns- The wave energy per length of field line is 
then proportional to 9^. Multiplying by the length ~ 3i?max of the excited field lines gives a 
net energy proportional to 9^ ~ -RNsZ-^max- 



The damping rate of this external torsional mode is examined in § [5 .21 , and compared 
with the observed peak luminosities of the giant flares. 



3.2.4- Constraints from the Rise Time 

The August 27th event showed a very steep rise in intensity, with width trise <^ 4 ms 
( [Hurley et al. 1999a , Mazets et al. 191?D ). If this timescale is limited by the propagation 
of a fracture, then a triggering zone moving with fracture speed Vf has size 

^ ~ Vfitrise = 4km(trise/4 ms). The fast rise of the 1979 March 5 event, trise^ 0.2 ms (Cline 
et al. 1980; or trise ~ 1 ms as reported by Fenimore et al. 1996) may require that a fast 
magnetic reconnection front was generated in the magnetosphere at the onset of the event 
(Paczyhski 1992). Alternatively, a fast rise could be understood if the large-scale, rotational 
motion results from self-organizing growth of collective motion in smaller-scale units in an 
"avalanche" effect. If these units have dimensions comparable to the depth of the crust, 
ARfj^ ~ 0.3 km, then rise times as short as ~ (Ai?^/V^) ~ 3 x 10~^ s are possible. 



3.3. Giant Flares in the Context of SGR Activity 



It has been argued that the giant flares involve a large-scale propagating fracture of 
the neutron star crust driven by magnetic stresses in the core (TD93; TD95). We have 
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described a specific configuration of the magnetic field, involving a strong twist in the core 
and crust, which can release ~ 10^^ ergs in individual increments as long as the yield strain 
in the gravitationally stratified crust is ipcr ~ 10~^. Transport effects such as ambipolar 
diffusion acting on a timescale comparable to the source age of ~ 10^ yr allow consecutive 
departures from magnetostatic equilibrium and repeated flaring activity, because unwinding 
of the core fleld is strongly inhibited by the compositional statiflcation of the neutron star 
interior. This is an efficient way to take magnetic free energy from its principal reservoir - a 
magnetic fleld anchored in the deep crust and core of a magnetar (see eq. [^), and convert 
it to observable photon emissions, while evading the tremendous neutrino losses which are 
inevitable for all dissipation occuring locally in the interior. 

The development of strong, localized shear in the exterior magnetic fleld allows the 
rapid dissipation of a large amount of magnetic energy outside the star, enough to power 
the observed flare emissions (eq. [|26|). In addition, a signiflcant fraction of the released 
energy may be deposited in a torsional Alfven wave in the core, which transfers energy only 
gradually through the crustal lattice to the exterior. Moreover, a predominantly rotational 
deformation of the crust will leave behind non-potential terms in the exterior magnetic 
fleld which support persistent electrical currents which dissipate in part through Compton 
drag off the ambient X-ray flux (Thompson et al. 2000). The dramatic change in the 
non-thermal persistent flux and pulse shape of SGR 1900-1-14 following the August 27 giant 
flare (Woods et al. 2000) provides direct evidence for a strengthening of these non-potential 
components during the flare. More generally, such a corona has been hypothesized as the 
source of the persistent non-thermal emission in the SGRs and AXPs (TD96; Thompson 
et al. 2001), and for the transient non-thermal emission following the 29 August 1998 
aftershock from SGR 1900+14 (§5.2.3 in Ibrahim et al. 2000). 

Shorter ~ 1 — 10 s bursts were observed soon after the two giant flares (see Mazets et 
al. 1979a for SGR 0525-66, and Ibrahim et al. 2000 for SGR 1900+14), with luminosities 
comparable to the pulsating tails of the flares. In the above model, these shorter bursts 
probably involve shorter fault-line slippage than the flares. The relation of the much more 
common ~ 0.1 s SGR bursts to the giant flares is perhaps more mysterious; while these 
events could also involve a propagating fracture, they may alternatively be driven by a 
more localized and plastic deformation of the crust (TD95). 

4. Soft Pulsations: the Light Curve 

The hard spike and the extended soft tail of the August 27 event released comparable 
energies, E ~ (0.5 — 1) x 10^*^ ergs, even though their durations differed by a factor of 10^. 
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This suggests that a large fraction of the outburst energy was injected during the initial 
~ 0.5 s. Whereas the enormous luminosity of the spike requires a relativistic outflow from 
the neutron star, the spectral stability and large amplitude oscillations observed during the 
last ~ 300 s of the outburst require instead that the radiated energy was confined close 
to the surface of the neutron star. This energy, if released into the near magnetosphere, 
forms a dense, thermal plasma of pairs and blackbody radiation with a temperature 
T = 0.72 {E/IQ^^ ergs)^/^ (L/10 km)-^/^ MeV (in a volume L^). The scattering opacity of 
electrons and positrons was then high enough to lock this energy onto magnetic field lines 
close to the neutron star, in a "trapped fireball" . This fireball cools by radiative diffusion 
through a thin surface layer, which contracts inward while remaining nearly congruent with 
the magnetic field (TD95). 

Strong evidence for this trapped fireball model comes from the faster-than-exponential 
decline in the X-ray flux at the end of the August 27 outburst, which is consistent with the 
hypothesis that the energy resevoir is evaporating completely in a finite time. Except for 
the first ~ 40 s (one tenth of the burst duration), the envelope of the August 27 light curve 
is well fit by the following function (Paper I) 



If the cooling luminosity is assumed to vary as a power of the remaining fireball energy, 
Lx oc E"", then x = a/(l — a). In the simplest case of a fireball with uniform energy density 
and surface energy flux, the index x is effectively the number of curved directions of the 
fireball surface: x = Dc = 2 for a spherical fireball; x = -De = 1 for a fireball filling bundle 
of flux lines that is (locally) cylindrical; and x = -De = for the case of a bundle of flux 
lines that is (locally) a thin slab (Fig. |^). 

This function (|3^ accurately captures the sudden final drop in flux seen in the Ulysses 
and BeppoSAX data (Paper I). The best-fit fireball index is x = 3 (or a = 0.75), and 
the evaporation time tevap = 375 s (Figs. 2-3 of Paper I). Note that the fit is excellent 
following the initial ~ 40-s of the outburst. In comparison, the best fit exponential profile 
exp(— t/texp) adequately describes the intermediate portion of the decay (Fig. 1. of Paper 
I), with a favored time constant texp = 78 s in the Ulysses data; but significantly overshoots 
the final ~ 40 s of the outburst. 

During the first ~ 40 s, the measured burst flux significantly overshoots both the 
exponential and trapped fireball fitting functions. The mild spin-modulation of the light 
curve during this smooth tail requires a more extended and variable photosphere, which we 
argued in Paper I is pair-dominated and results from a continuing creep of the neutron star 
crust in the active region (or possibly a standing torsional wave in its interior). The cooling 




(34) 
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of a continuously heated, pair-dominated corona is discussed in detail in §0. 



4.1. Fireball Index 

We now consider some simple physical models for a contracting fireball. Because the 
best-fit fireball index x is larger than 2 (the maximum number of curved directions for 
the fireball surface in three spatial dimensions), a homogeneous fireball cannot explain the 
data. In this section, we consider the effects of temperature and magnetic field gradients. 

The surface of the fireball will be assumed to have curved directions and 2 ~ Dc 
fiat directions, and to be fully symmetric with a single perpendicular coordinate r. Before 
the cooling wave propagates into the fireball, the temperature and magnetic field vary as 
powers of r: 

T oc ; 5 oc r^. (35) 

Temperature gradients perpendicular to B easily persist over the duration of the August 27 
event, because the photon diffusion time across a distance ~ i?NS through the fireball is 
much larger than the observed burst duration (or the time for the cool boundary layer to 
propagate into the fireball: TD95 §3.2 and §3.4). 

If the surface flux from a trapped fireball varies as 

FocT^'B^^. (36) 

then, because the emitting area scales as A oc r^'=, the luminosity is 

Lx = Lx{Q)[-) . (37) 

The energy density within the fireball includes contributions from both photons and 
electron-positron pairs. In two regimes of relevance, it has the form 

V o^T^'^B^\ (38) 

In particular, if many Landau levels are populated and T ^ me, then one recovers the usual 
expression for a relativistic pair plasma, V = (ll/4)aT^ where a is the Stefan-Boltzmann 
constant, and fii = 4 and fi2 = 0. The condition for a such 3-dimensional pair gas i^ 
2.7 T ;^ (2e5)i/2 equivalently 

/ ^ \ 1/2 

T ^ 0.3 MeV (39) 



5, 



QED, 



^In units with Ti = 1 = c, used frequently hereafter. 
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in a magnetic field B -Bqed (eq- [T|)- This condition is easily satisfied near marginal 
confinement, when the pressure of pairs and photons is close to the limiting value B'^/Sn 
(TD95). 

When the magnetic pressure greatly exceeds the thermal pressure, and T is not much 
less than rrie, the pairs dominate the thermal energy density: U ~ f/e+e- ~ (l/12)e-BT^ (in 
units where h = c = 1). The indices /ij are then /ii = 2 and /i2 = 1 [see eq. (53) and the 
discussion following eq. (56) in TD95; see also Figure 5 in Duncan 2000 and Kudari 1997]. 

The surface X-ray luminosity is related to surface area A of the fireball and its rate of 
contraction through Lx = —AU{dr/dt), or 

Equations (|37D and ( ^OD together give a simple differential equation for r{t). The solution 
implies that Lx{t) has the form of eq. (0) with 

l+7(/il-fTl)+/?(/i2-^2) ^ ' 

The case of a structureless trapped fireball, 7 = /3 = 0, reduces to x = -De < 2. Since 
X ~ 3 empirically, we infer that the fireball must have had structure, with indices which 
satisfiy 

(3/xi - 4ai)7 + (3/^2 - 4^2)/? = D,-3. (42) 

Specific values of the structure indices 7 and P, for several plausible assumptions about 
the fireball geometry and physical conditions, are given in Table 1. 

These values of 7 and (3 depend on specific choices for the parameters cTj in the surface 
X-ray flux (^) from the trapped fireball. One generally expects a positive dependence on 
temperature and magnetic field because of the thermal nature of the emission {ai > 0), and 
because Compton scattering is suppressed in a strong magnetic field ((T2 > 0). We have 
considered three cases which span the possible range of cooling-wave behaviors: one case in 
which F is most sensitive to B; one in which F is sensitive to T, and one with sensitivity 
to both variables. The case Ui = 0, (72 = 4/3 corresponds to a simple Eddington-limited 
flux from an atmosphere that is held down by gravity, when the magnetic suppression of 
scattering is taken into account (eqs. [45] and [93] in TD95). The case (Xi = 5/2, cr2 = 
corresponds to a magnetically-confined fireball in which the confined plasma fills many 
Landau levels (/ii = 4, /i2 = 0), and in which the integrated surface flux across B is limited 
by the rate at which energy is advected along the magnetic field, so as to smooth out 
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pressure imbalances (eq. [89] of TD95). Finally, the case ai = 3/2, (T2 = I (eq. [86] of TD95) 
corresponds to a similar process of advection-limited cooling, but in the opposite regime 
where the confined plasma is too cool to fill higher Landau levels (so that the radiative 
opacity across B is suppressed by the strong magnetic field). In these last two cases, T 
is interpreted as the temperature just inside the geometrically-thin radiative layer at the 
surface of the fireball. 

Table 1 shows that, in all physical regimes considered, a fireball that is slightly hotter 
in the center than near the edges, 7 < 0, can fit the observed light curve. Such a pattern of 
energy-loading in the magnetosphere could plausibly occur. An alternative, less attractive 
possibility, is that the temperature inside the fireball is uniform while the magnetic field 
declines shghtly toward its center, /3 > 0. 



4.2. Effects of Neutrino Cooling 

The trapped fireball light curve will also be perturbed by neutrino pair emission, 
e'^ + e~ ^ u + jy, at a rate U{e~^e~ vv) = 1.3 x 10^^ (T/MeV)^ ergs cm~^ s~^ (Dicus 
1972; Schinder et al. 1987). This can reduce the fireball energy on a time scale comparable 
to the duration of the giant outburst (§5 in TD95): 

- 28 ( — ^ ) s. (43) 



U{e+e-^uu) VlMeV 

This expression assumes that the temperature is high enough to excite pairs in the upper 
Landau levels (eq. ||3l 



The effects of pair neutrino cooling on a fireball light curve are shown in Figs. Hal and 



|B for several initial values of the dimensionless ratio (E/LxTu) and a fixed evaporation 



time tevap- In Fig- 0, we plot fireball light curves for various geometries {Dc = 1,2), 
assuming that the trapped pair gas is 3-dimensional {fii = 4, /i2 = 0) and that the cooling 
wave is limited by advection (ai = 5/2, a2 = 0; see eq. [89] in TD95). The neutrino 
cooling time is parameterized by the initial value of the dimensionless ratio E/ LxTy, where 
E is the trapped fireball energy. One sees that neutrino cooling introduces additional 
curvature in the (logarithmic) light curve, predominantly during the initial decline. This 
curvature is absent for a homogeneous fireball (7 = /? = 0), as well as for a fireball with 
powerlaw temperature and magnetic field profiles (constant 7, (3). In Fig. |?5|, we consider a 



combination of neutrino cooling and temperature gradients, with the value of 7 perturbed 
slightly from the best-fit value (^2|). (Here the initial value of {E/ LxTy) is calculated using 
the initial temperature at the outer boundary of the fireball.) 



-27- 



To summarize: neutrino cooling introduces curvature in the trapped-fireball light 
curve, which cannot be compensated by a power-law temperature gradient. A large value of 
[E/LxTy) ^ 10 does not fit the observed light curve, but one can not rule out a combination 
of slower neutrino cooling and more complicated temperature profiles. 



5. Bound on the Stellar Dipole Moment 

In this section, we quantify how the long duration and enormous peak luminosity of 
a giant flare set lower bounds to the magnetic moment of its source. First, by including 
the effects of neutrino cooling on the declining X-ray flux of a trapped pair-photon plasma, 
we are able to set a lower bound on Bdipoie that is less dependent on assumptions about 
the configuration of the field than those made in TD95. Second, by considering in more 
detail the coupling between a low-frequency internal (e.g. torsional) mode of the star and 
its magnetosphere, we are able to show that such a low-frequency mode cannot easily 
lose energy at the rate observed in the hard initial spike: even under the most optimistic 
assumptions about the damping of the external magnetic shear, the surface flux density 
must be very high (in excess of 10^^ G). 



5.1. Confinement of the Radiating Plasma During the Soft Pulsations 

A strong lower bound on the surface magnetic field is deduced from confinement of 
a large amount of energy, E ~ 10'^^ ergs, close to the SGR source during a giant flare. 
Assuming that the field is a centered dipole, and requiring that ^VStt > + Pe± at the 
outer bondary of the fireball, one deduces 



G (44) 



for a relativistic plasma (TD95). Here AR is the characteristic size of the fireball. 

The question which next arises, is how this bound is modified if the magnetic field 
is not a centered dipole. Indeed, the remarkable four-peaked pattern of the August 27 
pulsating tail provides direct evidence for higher multipoles (Paper I; Section H). Consider 
first how this bound is modified as the footpoints of the confining magnetic field are brought 
closer together. The minimum pressure of the confining magnetic field increases with 
decreasing footpoint separation L, S^j^/Stt ~ E/3L^ but the net magnetic moment of the 
confining field decreases^ jj^nin ~ -Bmin-^^ oc E^I'^E^I'^ . At the same time, the temperature of 
the confined plasma (assumed to have a 3-D distribution) increases, T oc 
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The bound on -Bdipoie cannot, however, be reduced arbitrarily by shrinking the firebalL 
The reason is that the neutrino coohng time of the fireball decreases rapidly with decreasing 

oc T-^ oc L^^/^ (45) 
at fixed energy E, even while the photon cooling time (eq. [Q) hardly changes: 

tevap ~ L/V,ooX ~ (46) 

in the absence of neutrino cooling. (The propagation speed of the cool boundary is related 
to the surface X-ray flux by Kooi = Fx{T)/[^aT^] oc T-^'^.) 

To get a sense of how far -Bdipoie can be reduced and remain consistent with magnetic 
confinement of the radiating plasma, let us consider a small offset dipole, centered at the 
outer edge of the star (Fig. Neutrino cooling has a large effect on the fireball lightcurve 
when Ty <^ tcvap, where the evaporation time can be written as 

tevap = (X + 1) (^) • (47) 

Indeed, in the case of rapid neutrino cooling, one can divide the cooling process into two 
phases: an initial neutrino-dominated phase where the fireball volume hardly decreases; 
and a subsequent photon-dominated cooling phase. Integrating the neutrino losses up to a 
time ~ tevap/ix + 1); the fireball pressure is reduced by a factor 



(Pe± + ^7)0 VT( 



4 



5 f E 



4 KLxTi,) Q 



-4/5 



(48) 



Thus, as E jLxTv pushed above unity, the lower bound on the magnetic moment 

does not change significantly: /^min oc [E^L^')^!'^ is almost independent of L. For the 
purposes of estimating /imin, we conservatively take {E / LxTy)^ = 10 in what follows. This 
corresponds to Ti, ^ 0.025tevap when x = 3. 

This lower bound on the neutrino cooling time yields an upper bound on the fireball 
temperature through eq. (^3]), 

In order to contain the energy radiated in the August 27 pulsating tail, the volume of the 
fireball must exceed 
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The confining magnetic field at its weakest point must be stronger than 
-^min/S^ = P-y + Pe± = if'^^^) which Combined with eq. (^9|) implies 
Brnin = 8.3 X 10^^{T/T^^^f (tevap/400 s)'^/^ G. Given that the plasma fills half the 
dipole magnetic field, and only along field lines which never drop below this flux density, 
the magnetic moment can be calculated to be 

327r y 10"^** ergs/ VTmax/ V400 s/ 

In this expression, -Bmin is the minimum flux density at the magnetic equator. 

This expression provides only a lower bound to the net dipole moment of the neutron 
star, for two reasons. First, no change was observed in the spindown rate of SGR 1900+14 
in the few months following the August 27 event (Woods et al. 1999), which implies that the 
dipole moment of the confining magnetic field was only a fraction of the total (Thompson 
et al. 2000). A conservative upper bound of ~ 30 percent to the fractional change in P 
translates into a maximum change 5/inet//^net ~ 0.15 in the magnetic moment in the case of 
magnetic dipole spindown (P oc /inet)- The tolerable change in the magnetic moment could 
be larger, (5/inet//^net ~ 0.3, if the spindown were accelerated by a persistent wind of particles 
and Alfven waves, assuming that the wind's power L\y does not change across the flare 

1 /2 

(P OC yUnet-Z^vK ! Thompsou & Blaes 1998; Harding et al. 1999; Thompson et al. 2000). On 
the other hand, the giant outburst is expected to produce a change 5/inet//^net ~ yUmin//^net 
in the net magnetic moment, and so we arrive at 

^ (S^MnetY^ O in32 /^'^/^net//inct E \ f T 3 

''"'^[^) "-"=2x1° i^^j [w^sjif^zJ 

(52) 

In this case, the effective polar dipole field is 

— tf^--"(lo4;s)"(£)"^(l^)"' O, ,53) 

A single dipole of magnitude (^2|) barely contains enough energy to power a single giant 
outburst, which suggests that either the entire external field of SGR 1900+14 was strongly 
sheared during the flare; or the actual field is much stronger than (|53|); or the surface of 
the star is covered with many dipoles ^ /Xmin- If the net energy output in giant flares is 
comparable to the external magnetic energy, then we estimate 

(54) 

(Note that this expression describes both a single large dipole, and also the incoherent 
superposition of A^ mini-dipoles of magnitude /imin; see also §12.3 in TD93.) 
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5.2. Peak Luminosity and Bound on the Surface Magnetic Field 

The initial spike of the August 27 flare reached a peak luminosity^ L ^ 10^'* ergs s~^ 
(Mazets et al. 1999), approaching a million times the Eddington luminosity of a neutron 
star. Such a high rate of release of energy allows us to set a stringent lower bound on the 
magnetic fleld in the source, which is compatible with the value deduced from the total 
expected flaring output (§^. Note that the opacity of the emitting material is not directly 
relevant to this bound, since the hard spike appears to have involved an expanding flreball. 

Motions of a neutron star crust will release magnetic energy directly by inducing 
tangential discontinuities in the fleld, which lead to fast reconnection and high frequency 
wave motion (Thompson 2000; §^.2.1| above). The output induced by a horizontal motion 



through a distance ~ -Rns over a timescale Atspike is i^spike ~ ex(-B^/87r)i?^g/Atspike. Here 
ex is the fraction of the magnetic energy which is dissipated and converted to high energy 
photons, and is estimated in eq. (E^). The implied r.m.s magnetic fleld is 



0.1 y V 10^^ ergs 



B^2x 10^^ — "-P'f ^P'"" G. (55) 



Much of the dissipated energy would initially be conflned by the (dipole) magnetic 
fleld. Nonetheless, the conflned pair plasma generated in such a giant flare can easily exceed 
the dipole pressure (which decreases as R~^) at a short distance above the surface of the 
neutron star. This provides a mechanism for driving a relativistic, pair-loaded flreball. The 
rise time of the August 27 hard spike could, in fact, be limited by the breakout of this 
conflned pair plasma. Evidence for repeated breakouts of pair plasma is provided by the 
~ 20 ms modulation of the hard spike of the March 5 event (Barat et al. 1983). 

Let us also examine the alternative possibility that the magnetic fleld merely acts as a 
couple between an internal mode of the neutron star and external Alfven-like excitations. 
Of the possible internal modes of a neutron star, a crustal shear wave has the strongest 
coupling to the magnetosphere (Blaes et al. 1989). The required magnetic fleld turns 
out to be similar, even after making the most optimistic assumptions about the damping 
rate of the external mode. This is because the fleld lines are anchored in very dense 
material, within which a shear disturbance propagates much more slowly than it does in 
the magnetosphere. (Equivalently, the equilibrium energy stored in the external mode is a 
tiny fraction of the energy in the internal mode; § ^.2.2| .) A related mechanism involving the 



^This applies after allowing for dead-time corrections; this value remains below that inferred previously 
for the peak luminosity of the March 5 event at the distance of the LMC (Mazets et al. 1979). 
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coupling of an internal p- or /-mode to the magnetosphere was examined soon after the 5 
March 1979 flare, by Ramaty et al. (1980) and Lindblom & Detweiler (1983). As we argue 
below, the damping rate was significantly overestimated by these authors, and is in fact 
much closer to the rate for vacuum dipole radiation. 

Global toroidal modes in the crust can be excited directly through a large-scale 
fracture, as is expected in a giant fiare, or indirectly through coupling to a torsional Alfven 
mode in the core (§ |3.2.2| ; Duncan 1998). The fundamental toroidal mode (with n = 
radial nodes) has a frequency u ~ + l)(y^o/ Rt<ss) for ^ > 2 (McDermott et al. 1988). 
Whereas the torsional Alfven mode has a frequency u = nVA^z/'^Rf^s (in a uniform poloidal 
magnetic field with corresponding Alfven speed Va,z = -B^/V^vrp; eq. [0). To begin, 
let us consider an excitation confined to the crust. Suppose that an energy E'eias ~ 10'^'' 
ergs was deposited in toroidal modes (of various harmonics i) during the initial ~ 0.5 s 
spike. The r.m.s. mode amplitude (^^)^^^ at the neutron star surface is releated to the 
mode energy i^eias by substituting ip^ ^ i{i + 1)(^/-Rns)^ in equation (H), which gives 
^eias = liii + 1) McouiombV^^o / Rf,sf . Wc havc in turn 

f^^^^^H^V [ ^ km. (56) 

VO.O2M0/ VlO^cms-V ^ ' 



2\l/2 



0.1 



elas 



10^4 ergs 



+ 1) 



The resulting harmonic displacement of the neutron star magnetosphere can be divided 
into two components (§ |3.2.3|) : a quasi-static deformation of field lines which close within a 
radius -Rmax ~ |c(27r/co'); and an Alfven mode which propagates freely along more extended 
field lines, with a group velocity (icj/rfk = cB. 

To set a conservative lower bound on the rate of dissipation of energy in the 
magnetosphere, let us assume that the quasi-static twist is damped within a substantial 
solid angle Afi on a short timescale 1t[ juj '^^ tdamp ^ -Rns/c (Notice that the light-crossing 
time -Rns/c is much shorter than the mode period, by a factor ~ V^o/c-) The amplitude 
of the deformed field grows at a rate ddB-f^a^/dt ~ (ci;^/i?Ns)-BNS) leading to an equilibrium 
amplitude (5-Bns/-Bns ~ (ctdamp/-RNs)(i^^/c). The torsional mode energy then decreases at 
the rate 

_^.,^.«.,)|^e(^)^(f|=). (57) 

This expression is maximized at a damping time tdamp ~ -Pmode = 27r/u; (above which 
the damping rate scales inversely with tdamp)- It can be inverted to yield the minimum 
mode energy needed to power a given rate of dissipation, which we set equal to the spike 
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luminosity L^pii^e- 

r-, 1 in43 -^spike \ / -Pmodc \ f B^S \ ( ^Coulomb \ f A^7^ ^ _^ 

^-- = ^^^^ 1, 1044 ergs s-0 ll^J iomSMejtj ^^^^ ^ • 

(58) 

One sees that even if the mode is damped over a large solid angle AVt ~ 4-71, a surface 
magnetic field of at least 5ns 

~ 3 X 10^3 (Pmode/0.03 s)-i/2 Q jg required. Otherwise, the 
mode must carry much more energy than is radiated. Indeed, crustal shear oscillations are 
barely capable of storing sufficient energy (eq. |lTl|). 

It is perhaps more natural for a large energy to be deposited in a torsional Alfven 
mode in the liquid interior, for the reasons discussed in § |3.2.2| . The mode energy would 
then be even larger, -^Alfven ~ \M^q,{^uj pawen)"^ ■ In order to power the hard spike via Alfven 
waves from the excited star, a condition similar to eq. (|58D must be satisfied, but with the 
mode mass Mcouiomb replaced by Mns- The lower bound on the magnetic field is then more 
severe, 5ns 3 x 10^^ (P„,ode/0.03 s)-V2 g. 

These bounds on the mode energy are conservative: it is less than clear that a smooth 
periodic deformation of the external magnetic field will damp in a time as short as Pmode, 
and only a small fraction ~ 3i?Ns/c-Pmode of the external field lines are able to support a 
propagating shear Alfven mode with a frequency comparable to that of the internal mode. 
By contrast, dissipation will occur effectively outside the star through the build-up of 
strong, localized shear (e.g. near the site of a fracture). 

Could an internal toroidal mode power the more modest output of Lx ~ 10^^ ergs 
s~^, above and beyond the trapped- fireball light curve, which is observed during the 
initial ~ 40 s of the August 27 fiare? We ascribed this excess emission to a continuously 
heated pair-dominated corona in Paper I. (Fatuzzo & Melia (1993) and Melia & Fatuzzo 
(1995) have explored the related possibility that short SGR bursts may be powered by 
the continuous dissipation of sheared Alfven waves propagating away from a more weakly 
magnetized neutron star, but did not address the relation between the Alfven wave fiux and 
the mode energy.) Setting the mode energy equal to this excess X-ray output of ~ 10"^^ ergs, 
one deduces from equation (^) that the external Alfven mode must be damped through 
a solid angle Af] ~ 47r (5ns/10^^ G)'^ (p^^^^^/o.03 s)-^ (M^ode/MNs) Sr. For the reasons 
just described, this angle is uncomfortably but not impossibly large. The excess emission 
during the first ~ 40 s could, alternatively, have been powered by a continued shearing of 
the external magnetic field, driven by a diminishing creep of the crust in the active region; 
or by the gradual decay of static magnetospheric currents. 

Finally, let us compare our results with the model of Ramaty et al. (1980) and 
Lindblom & Detweiler (1983). It was suggested in these papers that an internal /- or 
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p-mode would be damped by shock-heating of the surface layers of the neutron star, and also 
by a strong flux of Alfven waves into the magnetosphere. Repeated shocking of the surface 
layers of a star can occur only if the time for a sound wave to cross a pressure scale-height is 
longer at shallower depths. (Otherwise, an upward-propagating sound wave will be reflected 
downward before it has a chance to steepen into a shock.) This condition is indeed satisfled 
in the envelope of a pulsating AGB star, but it fails to be satisfied in the crust of a neutron 
star. We start with the equation of hydrostatic equilibrium, dP/dr — —pGM[< r)/r'^, 
where P{z) oc [p{z)Y at depth z = i?NS - r. (The index 7 ~ | for p ~ 10^ - 4 x 10^^ g 
cm~^ and 7 ~ | at lower densities.) Integrating gives P{z)/p{z) ~ (7 — 1) 2; (GMnsZ-Rns)) 
along with the sound speed Cg ~ ^J^P/p and the pressure scale-height ip — P/ {—dP/dr) as 
functions of depth. In a neutron star, the sound-crossing time decreases outward. 



and so the surface layers are manifestly stable to repeated shocking. 

Now let us consider the electromagnetic damping of an /- or p-mode. The damping 
rate into Alfven waves is even lower, because an /- or p-mode does not shear the external 
magnetic field lines. The mode could have a frequency as high as a; ~ 10^ Hz, and so 
a much wider bundle of external field lines would be excited than by a lower-frequency 
torsional mode. (We assume for now that the magnetosphere is able to respond to the 
surface motions as a conducting fiuid.) To set a conservative upper bound on the energy 
fiux into the magnetosphere, we allow all the field lines to carry away Alfven waves of an 
amplitude 55ns/-Bns = C^/c and a surface energy flux c((5i?Ns)^/87r. The mode energy 
-E'mode — |-^NS<^^(^^) IS then damped at a rate 

_ ti-gmode ^ , d2 -^NS f ^^(^^) ] ^ -^mode-RNS-^NS /gQ\ 

dt ~ Sn \ J Mnsc ' ^ ^ 

For a neutron star mass M^s = 1.4 and radius i?NS — 10 km, this expression becomes 

dE^ode , .ri42 ( -^mode \ / -^NS ^ -1 

; ~ 1 X 10 — — ; ^ ergs s . (61) 

dt \1044 ergsy VlO^^ GJ ^ ^ ^ 

As a check, consider the damping rate due to vacuum dipole radiation. As before, we 
normalize the mode energy £^mode — I-^ns^^*^^ to the observed X-ray output of the flare. 
One has ^ 

_ f dE^node\ _ 2a/ 2 [ }_ f :^mode\ -^nS-^NS^^ /^on 

1, dt Jmdr"3^^ UnsJ -3 1mnsJ c3 ' ^^'^ 
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where M. — I^ns-^ns magnetic moment and B^s the polar surface magnetic field. 

For a mode frequency u; ~ 10^ Hz, M^s = IAMq and it^NS — 10 km, this expression 
becomes 

Using either calculation, the rate of dissipation fails by orders of magnitude to accomodate 
the observed output of the spike - even if the surface field is in the magnetar range. 

6. Flcire Spectra and Spectral Evolution 

6.1. Spectrum of the Hard Initial Spike: An Expanding Pair Fireball 

The first part of the August 27 event had a much harder spectrum than the remainder 
of the outburst. A very hard component dN/dE oc was measured by BeppoSAX 

during the first 67 seconds, with a much softer power-law component dN/dE oc 
remaining during the following 128 seconds (Fcroci ct al. 1999; Paper I). The detection by 
Konus/Wind of a dramatic drop in the > 250 kcV emission following ~ 0.7 sec post-trigger 
(Mazets et al. 1999) suggests that the very hard spectrum component was mostly emitted 
during the initial spike. 

The large-amplitude fluctuations detected in > 250 kcV photons at the end of the hard 
spike are also remarkable (Mazets et al. 1999). At 0.2 — 0.7 sec after the event's onset, 
these fluctuations were were more profound in > 250 keV photons than in the 15-250 keV 
band, and occured on timescales as short as ~ 10^^ s. This behavior is reminiscent of 
high-redshift, classical gamma-ray bursts (GRBs), in which the light curve is often much 
smoother at lower energies (Pendleton et al. 1997). 

The peak luminosity of the August 27 spike (L > 4 x 10^^ ergs s~^; Mazets et al. 1999) 
is intermediate, on a logarithmic scale, between a thermonuclear flash and a high-redshift 
GRB. But we have much more rehable, direct information about its source, SGR 1900-1-14, 
than we do about the sources of cosmological GRBs. Thus, detailed studies of the August 
27 flare may pay dividends in understanding the more mysterious classical GRB sources, 
even though the integrated spectrum of the August 27 flare is not typical of a classical 
GRB. 

The volume of the emission region is strongly constrained by the rapid variations in 
the > 250 keV flux. Although the hard power-law component is detected only up to ~ 700 
keV (Feroci et al. 1999), there is no evidence for a high energy cutoff. Therefore, let us 
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consider the consequences of a source function that extends above an energy ~ mgC^ in 
the (relativistically) expanding frame of the outflowing particles and magnetic field that 
power the hard spike. Pair creation will attenuate the flux above a measured energy FrngC^, 
where F is the bulk Lorentz factor of the outflow. Because the photon flux is dominated 
by high-energy photons {dU/dE oc E-^l^), one expects copious pair creation, and the 
luminosity in relativistic pairs is comparable to the photon luminosity. Pairs then make the 
outflow very thick to scattering, out to a considerable radius 

(64) 



\47rrmeC^ 

Here, Ai? ~ LAt is the energy carried by an individual spike of peak luminosity L and 
width At. In order for the spike not to be smeared out as the hot shell of ejecta expands to 
a thickness ~ i?/2r^, one requires cAt > i?T-=i/2r^, or equivalently 

r>r..in = 7 ^ \ (-^) . (65) 

ylO*'^ ergs s ^ ) VO.Ol s/ 

Most of the dissipation that produces the hard non-thermal spectrum must therefore occur 
at a large radius, 

fi>i?.=i^3xl0^° f \ (-^) (J-) cm. (66) 



10^3 ergs s-iy VO.Ol s/ VP 



mm 



The emission region of the non-thermal spectral component is very large, but 
comparable to the speed-of-light cylinder of SGR 1900+14 {cP/2ti = 2.5 x 10^° cm). It 
sits far enough out that direct synchrotron cooling off the dipole field of the source can be 
neglected — even if the field lines are combed out into a ~ geometry by the pressure 
of the escaping particles and radiation. By contrast, the ejecta from the outburst could 
carry a signficant magnetic field. For example, B ~ if F grows linearly with radius - as 
is expected in an expanding fireball which contains a magnetic field but whose pressure is 
predominantly thermal at the source. Alternatively, inverse-Compton cooling by advected 
X-ray photons will still be a powerful coolant at a radius ( pB[ ) . The cooling time of a particle 
of random energy '-^jrieC^ in the bulk frame is tcompton ~ 7~^(crTL/47r_R^F^meC^)~^, and so 
the ratio of the cooling time to the flow time 

Won . /_^\ (^\"'\ (67) 



R/cV \Rr=lJ VP 



mm 
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6.2. The Pulsating Tail Spectrum: Temperature Regulation through Photon 



Wc now consider the temperature of the radiation escaping from a trapped fireball, and 
compare it with the best (blackbody) temperature fit to the pulsating tail of the August 27 
flare. The temperature appears to have remained remarkably constant during the phase 
of large- amphtude pulsations, even as the X-ray flux continued to decline (Mazets et al. 
1999). The best two-component spectral flt (black body -|- powerlaw) yields kT^b — 10.8 
keV during interval B (Table 2 of Paper I) and 12.2 keV during interval C (Table 3 of Paper 
I). These values are consistent with the minimum photospheric temperature of a trapped 
fireball in super-QED magnetic fields, as calculated by TD95. 

Near the photosphere, the spectral shape is determined by two coupled processes: 
Compton scattering and the creation of new photons through splitting 7 ^7 + 7 (TD95). 
The photons flowing out from the base of the trapped fireball are expected to dominate the 
specific heat and so, in the absence of an extended pair corona, the mean energy per photon 
will remain approximately constant outside the fireball photosphere. Within the outer 
layers of the optically thick fireball, the photons can maintain a Planckian distribution at 
temperatures well below the internal fireball temperature of ~ 1 MeV. The mean energy of 
the escaping photons depends directly on the splitting rate as a function of frequency. In 
marked contrast with the strong scaling of the splitting rate in sub-QED magnetic fields, 
the splitting rate approaches a S- independent value in fields much stronger than Bqeb, 



(Adler 1971; Thompson & Duncan 1992). This implies immediately that an E-mode 
photon propagating a distance i?Ns ~ 10 km through a super-QED B-field will split if 
huj > 38(i?Ns/10 km)^i/5 j^^y (TD95; Baring 1995). 

Compton scattering becomes strongly anisotropic in a background magnetic field, 
with a frequency- dependent cross-section. Near the surface of the star, the energy of the 
first Landau excitation [about (2S/Sqed)^^^ '^eC^] is much higher than the temperature 
of the emerging X-rays. In this situation, there is a strong suppression of the E-mode's 
scattering cross-section: aE — {comec/eBo)^ ax (e.g. Herold 1979). By contrast, the 0-mode 
scatters with a cross-section near Thomson (except for propagation almost parallel to the 
background field). This suppression of the E-mode scattering opacity greatly increases the 
radiative transport rate close to the neutron star. One has 



Splitting 




(68) 




(69) 
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when the dielectric properties of the medium are dominated by vacuum polarization 
(Silan'tev & lakovlev 1980; TD95). First considered as a way to enhance emission from a 
hydrostatic atmosphere that is confined by gravity at the surface of an ultra-magnetized 
neutron star (Paczyhski 1992; Ulmer 1994; Miller 1995), the suppression of 15- mode 
scattering has potentially more direct applications to radiative transport across the 
confining magnetic field lines of a trapped fireball (TD95). 

The two polarization modes are also distinguished by their ability to split. Only the 
13-mode can split because - when vacuum polarization dominates the dielectric properties 
of the medium - only the energy and momentum of E-mode photons can be conserved by 
dividing into two obliquely propagating daughter photons. Splitting occurs predominantly 
via 

E^O + 0, (70) 

with E ^ E + O having a lower rate (e.g. Berestetskii et al. 1980). When the photon 
distribution is close to Planckian, net transport out of the E-mode by splitting is suppressed 
by the inverse process of merging O + O ^ E (TD95). A cascade of photons from high to 
low X-ray energies cannot occur purely through splitting, as conjectured by Baring (1995), 
even in magnetic fields as strong as ~ 10^^ G. 

Even in the region where the E-mode is able to stream freely, the 0-mode can still 
undergo many Compton scatterings and relax close to a Bose-Einstein distribution. This 
permits a very simple generalization of the LTE diffusion formalism to an anisotropic, 
magnetized plasma at large E-mode scattering depth, in which the photon energy and 
number fluxes are expressed as linear superpositions of gradients in the temperature T and 
photon chemical potential /i (TD95). As a result, there is a critical temperature above 
which the distributions of the E- and 0-modes both become thermal, which works out to 

(eq. [133] in TD95). This value matches well the best fit temperatures for the August 27 
event. The mean energy per escaping photon cannot fall below ~ 2.7kTsp, in the absence of 
new photon creation processes or strong adiabatic cooling. 

Now let us consider the spectral variations within each 5.16-s rotation period. The 
emission is harder at the peaks than in the valleys (see the hardness ratio plots in Fig. 7a 
of [Mazets et al. 195^) . A positive correlation between hardness and intensity is expected 



of thermal emission from an optically thick plasma. Nonetheless, the large-amplitude 
oscillations show overall an anti-correlation between hardness and intensity (Paper I). A 
deep trough is present in each rotation cycle after ~ 40 s, and plausibly involves occulation 
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of the trapped fireball by the neutron star. The spectrum evolves from hard to soft across 
the trough, which contributes to an overall sawtooth-like pattern in the hardness ratio (see 
Fig. 1 in p^'eroci et al. 1999] ). This hard-to-soft evolution suggests a significant deviation 
from mirror symmetry in the confining magnetic field with respect to the center of the 
trough. 



7. The Smooth ~ 40 Tail: A Continuously Heated Pair Corona 

The giant flares share two main features: an extremely bright, hard, and brief ~ 0.5-s 
transient at the onset of the burst; and an extended train of softer, large- amplitude pulses 
which repeat coherently at the rotational frequency of the star. The initial spike has the 
properties expected of an expanding, pair-dominated fireball; and, as shown in Paper I, 
the declining amplitude of the pulsations in the August 27 flare can be well fit by the 
contracting surface of a trapped fireball. 

Nonetheless, this simple picture of a sudden release of thermal energy within the first 
~ 1 s, part of which escapes directly and the rest of which remains confined close to the 
star, does not appear to provide a complete description of the August 27 flare. Within the 
first ~ 40 seconds, the X-ray flux declines smoothly and shows only a mild modulation at 
the 5.16-s spin of the source (Figs. 6 and 7 of Paper I). The flux exceeds the favored trapped 
fireball model during this interval (Fig. 3 of Paper I), and the spectral hardness decreases 
by a factor ~ 1.5 before flattening out after ~ 40 s (Mazets et al. 1999). This behavior 
was ascribed in Paper I to an extended, pair-dominated corona that scatters and heats the 
thermal radiation emerging from the inner, trapped fireball (Fig. |^). 

This ~ 40 s coronal component of the August 27 flare emitted about 20 percent of the 
total energy (excluding the hard ~ 0.5 s initial spike; §6.3 of Paper I). At t ~ 15 s, the total 
luminosity was twice the envelope of the trapped fireball, which according to eq. (|34D at 
that time hardly deviates from its initial value -Lx(O) = (1 -|- x) -^fireball /4vap- Here x — 3 is 
the fireball index, -Efirebaii the initial fireball energy, and tevap — 375 s the time at which the 
fireball evaporates. 

The cooling time of this corona is very short given the high burst luminosity. 
For an optical depth to scattering Tea ;^ 1 (as required by the radiative model 
discussed below) the cooling time is set by radiative diffusion to a magnitude 
~ Tes-Rcorona/c = 10^^ (Tgs/lO) (-Rcorona/30 km) s. Thus, the coroual heating must be 
almost continuous. It could be powered in at least two ways: through a persistent creep 
of the crust at the heated fracture site, which results in continued shearing of the external 
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field; or possibly by a coupling of internal shear oscillations to the magnetosphere (§ p.2.2|) . 
The flare light curve does not provide a simple diagnostic between these alternatives. 

7.1. Steady Compton Heating of the O-mode 

Why should a continuing release of seismic energy lead to the formation of a distinct 
scattering atmosphere, instead of merely adding energy to a trapped fireball which forms 
during the first ~ 0.5 s of a giant flare? We now show that there is a critical rate for 
injection of energy into the trapped fireball, above which the energy becomes trapped before 
it can leak out radiatively; but below which a steady balance between heating and diffusive 
radiative cooling is possible. This critical luminosity depends on the nature of the photon 
source. When double-Compton emission dominates, the answer is given in Thompson 
(1997). In the present context, photon splitting 7 — > 7 + 7 is the most effective source of 
new photons above a critical temperature T ~ 11 keV (TD95) 

We first summarize the properties of a pair corona with a temperature kTe <C mgC^ in a 
very strong magnetic field B ^ -Bqed- A realistic heating mechanism involves electrostatic 
acceleration of the pairs through a turbulent cascade of interacting Alfven waves: the 
current density is driven just above the value that can be supported self-consistently by 
the available particles, and a displacement current is induced parallel to the background 
magnetic field (Thompson & Blaes 1998). 

The pairs cool primarily by Compton upscattering the O-mode, due to its much larger 
cross section. Although at large scattering depth radiation diffuses fastest via the El-mode, 
the orthogonal O-mode will be regenerated^ by scattering (with cross section teo — te'-, 
Meszaros 1992) and by splitting (TD95; Miller 1995). Moreover, in some circumstances 
the optical depth may be too small to allow effective mode transfer O ^ E hy scattering. 
(Conversion of the O-mode to the E-Yoode does not occur by splitting, except at enormous 
optical depths where the plasma dominates the dielectric properties of the medium.) 

In the following discussion, we will assume that the electrostatically heated pairs are 
sub-relativistic and pinned in the lowest Landau level. Indeed, at the very high compactness 
(£ = Lxc^T / '^T^Trif.c^ R^Q, ~ 10^) characteristic of the August 27 smooth tail, the pairs will 
equilibrate at a temperature close to that of the O-mode photons. Since the effective 



^Indeed, the inability of high energy O-mode photons to spht, combined with their rapid Compton 
heating by pairs, provides a mechanism for generating the non-thermal X-ray spectra of SGRs in their 
quiescent (non-bm'sting) state. 
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temperature of the photons is kTo ~ 10 — 20 keV, the spectral distribution of Comptonized 
photons will have an exponential cutoff F,^ oc exp(— /izz/fcTg) at high energies hv ^ kT^. 

The question which we address here is: what is the range of temperature for which 
the pair plasma is stable to an upward excursion in the density 11^+ + rig- ^ 211^+ of 
electrons and positrons? An accompanying increase in the pair temperature Te forces a 
runaway of the pair density (due to the higher flux of Comptonized photons above the 
pair-creation threshold hu ~ mgC^). This runaway is inevitable if the pair plasma is in local 
thermodynamic equilibrium (TD95). 

If, instead, decreases when ne+ increases, then a stable balance between electrostatic 
heating and diffusive photon cooling is possible. Then the distribution functions of E-mode 
and 0-mode are both approximately Wien. Such a steady balance requires a continuous 
source of 0-mode photons, which can be provided either internally by splitting the E-mode, 
or externally by scattering incident E-radiation to the O-mode.^ Splitting can be expected 
to dominate where the magnetic field is stronger than -Bqed! but an external source of 
E-mode photons is necessary in a more extended corona around a bursting magnetar. We 
consider each case in turn. 

7.2. Equilibrium Temperature of the Corona 

We first consider the equilibrium temperature of the corona, given a steady balance 
between the rate Lcorona at which the pairs are electrostatically heated, and radiative 
transport out of the magnetosphere. The result depends on the mechanism by which 
radiative transport occurs (through the E-mode or the 0-mode). We idealize the corona as 
a sphere of radius R. 

If the source of photons lies below the corona - i.e. in the form of a cooling trapped 
fireball - then photons diffuse into the corona predominantly via the E-mode, before being 
converted to the 0-mode and heated. One deduces that teo ^ 1 (and that energy is also 
transported out of the corona through the E-mode). In addition, the scattering depth to 
the 0-mode is very large, 

( ^ y 

Teo 57r2 \5qed/ \meC^ 



^Double Compton scattering of the O-mode — O + — > O + O + — is an effective source of new 
photons only at larger scattering depths than are needed to convert the E-mode. 
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(Here Te is the temperature of the external photon source, and we have made use of eq. 
|69| .) Fresh 0-mode photons are rapidly Compton heated up to a Wien distribution with 
temperature Tq — and mean energy per photon SkT^. Furthermore, the E-mode and 
0-mode photons are in equihbrium at large teo, and one has Te — T^, he — no- 



The balance between heating and diffusive radiative cooling becomes 

-^corona 3kTEnE 



AttR'^c 3{2ne+aE{TE)R) 

As just described, such a steady balance requires that photon creation be too slow to 
establish local thermodynamic equilibrium. The densities of photons and pairs are then 
related by 



(73) 



ne+ no 



ne+{LTE) no{LTE)' 
The corresponding LTE densities are 



(74) 



and 

2no{LTE) = no{LTE) + nE^LTE) = 0.244 ( ^ J (76) 

at temperature Tg. The equilibrium temperature is then determined from eqs. (|69D and 
(0), 

-^corona / \ / TflpC \ ( B \ f R 



1034 ergs s-i V^ecV ''''^ 1^ fcT^ J I^SqedJ UO kmj ' ^^^^ 

The temperature corresponding to the brightest, short SGR outbursts (~ 10^^ ergs s~^) in 
super-QED fields {B/Bqeb ~ 10) is /cTg ~ 25 keV, which increases only to ~ 50 keV at 

i^corona ~ 10=^^ CrgS S'^ (Fig. §) . 

It is also possible that the photon source is internal to the corona (e.g. splitting). The 
corona can then maintain a much lower scattering depth while cooling diffusively through 
the 0-mode - without any need for converting photons back to the E-mode. In this case, 
the balance between heating and cooling becomes 

-^corona SfcTo^o x 

AttR'^c ^ 3(2ne+aTi?)' 

since cxq ~ ctt after averaging over angles. The equilibrium temperature is now determined 
by 

kT, (m,c^\ ( B \ ^ / R 



1034 ergs s-i Uec^ ) ''''^ V ^T, ) Uqed j UO kmj ' ^^^^ 
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Notice the different dependence on both Te and B/Bqeb- In this case, the temperature 
corresponding to ^corona ~ 10''^ ergs s~^ in a magnetic field ~ 10 5qed is shghtly lower, 
kTe ~ 20 keV, increasing to kT^ ~ 30 keV at Lcomna ~ 10^® ergs s~^ (Fig- ^• 



7.3. Critical Luminosity: Photons Created by Splitting 

We now look for the maximum energy input to the corona that allows a steady balance 
between heating and diffusive cooling. The key criterion is that the equilibrium pair density 
ne+ be smaller than the LTE value i^^. 

Photon splitting will occur within the bulk of the corona. We can assume that photon 
transport out of the corona is through the 0-mode, because splitting converts the E-mode 
to the 0-mode (but not the reverse). As a result, transport cannot occur self-consistently 
via the E-mode. A high temperature kTo ^ 20 keV is needed to generate sufficient pairs 
to scatter the 0-mode back to the E-mode {tqe ^1); but because this temperature is 
well above the critical value (~ 11 keV) where splitting is rapid, the transport of E-mode 
photons will be limited by splitting, Lcorona/47r i?^c ~ SkTEUEc/TspR. In this situation, it 
would not be possible to maintain a steady balance between the creation and the outward 
diffusion of photons, because the rate of creation of photons ~ AnnER^Tsp would exceed 
Lcorona/^kTE by a factor ~ (Esp-R/c)^ > 1. 

Instead, the radiative flux out of the corona is carried by the 0-mode. A steady 
balance between photon creation and loss requires that each diffusing 0-mode photon 
have a substantial probability (~ |) of converting to the E-mode and splitting.]^ Thus 
teoTo ~ 1, which corresponds to 



QED, 



But 



ro^{2n^.)arR = l.^{-^^ (lO^^^i^) (2^) [w^] ^''^ 



(making use of eqs. |]7S[ and |]7U|). The bound ne+ < ne+{LTE) then implies, together with 



^"^Once again, we must assume that the photon density is well below the black body value value at 
temperature Tq to obtain a steady balance between heating and cooling. Thus, the inverse process of 
photon merging can be neglected. 
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eq. dg), 

Leorona < L™(max) = 4 X 10- (^) (^) ergs s^. (82) 

A stable balance between heating and diffusive cooling is possible below the luminosity 
-^corona (max). In such an equilibrium, the electron temperature is strongly buffered and 
remains close to 20 keV as Lcorona varies by a few orders of magnitude around ~ 10^^ ergs 
s~^ (eq. [0). The right side of eq. ( P^ scales as ~ i? in a dipolar magnetic field. 

This equilibrium between heating and cooling is steady, as it is straightforward to 
check. A rise in 71^+ implies an increase in tq, which must be compensated by a decrease in 
teo to maintain the relation tqTeo ~ 1 (eq. |80|). By contrast, when the fireball is in local 
thermodynamic equilibrium, an increase in ne+ is directly tied to an increase in T. In such 
a situation, continuing energy input will force a runaway increase in the energy density to a 
trapped fireball (TD95). 



7.4. Critical Luminosity: External Source of Photons 

Now consider the case where the corona is sufficiently extended that the magnetic 
field is weaker than -Bqed, and photon splitting can be neglected over most of its volume. 
Diffusive cooling now occurs through the E-mode, and a steady balance between heating 
and cooling fixes the electron temperature through expression (^). Given a source of 
fresh E-mode photons (of temperature TJ^) escaping from the trapped fireball, the optical 
depth through the corona is fixed. In the regime of interest, teo > 1? and the soft fireball 
photons are heated by only a fraction ~ t^q of the coronal input Lcorona: the equilibrium 
pair temperature is fixed by energy conservation. 

X 3(Ts - T|") = ^^£^. (83) 

(We neglect the difference between the mean energy per photon in the blackbody and Wien 
distributions. This expression is approximately valid even for low teo, where only a fraction 
~ Teo of the fireball luminosity is captured by the corona.) Combining this relation with 

-^^corona = {'ikTE)Nf givCS 

TEO ^ (l - ^) = (84) 

At the high luminosities of SGR outbursts, the pairs in the outer corona equilibrate 
to a temperature only slightly different from that of the heated 0-mode photons. During 
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repeated scatterings, the mean energy per photon exponentiates as e^, where the Compton 
parameter is 

y ~ — . (85) 



In equihbrium, 
which imphes 



l/ = ln(l + |^), (86) 



^fireball 



To-'-iE[^)'''- '''' 



We can now determine the critical luminosity below which a steady balance between 
heating and cooling is possible. One has, from equations (|T3D and (fTSp, 

teo ~ {2ne+)(JE{TE)R = 1.5 



» / -'^corona \ / '^-'■e » / -'^^ 



ne+{LTE) J VIO^^ ergs s-i; \20keVj VlO km 

(88) 

This expression is very similar to (pij), with the distinction that diffusive transport is 
through the E-mode instead of the 0-mode. As a result, Te is slightly higher than in the 
case where fresh photons are created by splitting (but is still strongly buffered through eq. 



|77|). The critical luminosity is determined by setting ne+ < ne+{LTE): 



42 /i '^E^\ ( \ f R \ 1 



i.„o.. < i„„„.(max) = 1.5 X 10- [I - j ^— j erg. s"'. (89) 

As before, this equilibrium is steady as long as the densities of 0-mode and (heated) 
E-mode photons lie below the blackbody value at temperature T^ — Tq. That is because a 
rise in ne+ increases teo^ and hence reduces the equilibrium (Wien) temperature to which 
the escaping fireball photons will be heated (eq. p3|). 



After the coronal heating rate ^corona drops below the luminosity of the cooling fireball, 
this pair atmosphere evaporates and the scattering photosphere contracts to the outer 
boundary of the trapped fireball, where the opacity is dominated by ion-electron plasma 
(TD95). This provides an explanation for the flattening of the light curve simultaneously 
with the appearance of large-amplitude oscillations. We outline the effects of the strong 
magnetic field on this Compton corona in the next section. 
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7.5. Very High Current Densities 

In the preceding, we have assumed that the pairs are heated electrostatically by a 
fluctuating current, and derived the critical heating rate above which the pair-photon 
plasma runs away to local thermodynamic equilibrium. There is, similarly, a critical 
current density J above which photons will remain trapped in the current-carrying region. 
The energy density of the photons exponentiates, thereby creating a pair plasma in local 
thermodynamic equilibrium. 

The optical depth to the 0-mode through the current-carrying charges is 

To ~ {ne- + ne+)crTR ~ ^(Tr-R, (90) 

where R is the size of the current-carrying region. In a pure pair plasma, with equal 
numbers of charges moving in opposite directions along the magnetic field lines, the mean 
frequency shift per scatter is second order in the drift speed V of the current carriers, 
Au/u ~ [V/cy. In equilibrium, the temperature of the photons and (1-dimensional) pairs 



IS 



A;To^(^)'. (91) 

When To > I, each 0-mode photon undergoes ~ Tq scatterings, and the Compton 
parameter 

depends only on the current density J and the size R of the 'corona'. However, Tq is 
strongly buffered by pair creation and annihilation (§ [7.21 ). Thus, above a critical current 
density J, the scattering depth becomes too large to allow the photons to diffuse out before 
their number and energy density multiplies. For example, in the case where new photons 
are created by splitting, the critical current density corresponds to an optical depth tq 
equal to (^. 



7.6. Alternative Models for the Smooth Tail 

Let us consider two alternative explanations for the smooth tail: a steady increase 
in the ion-electron loading at the fireball surface; and neutrino cooling. Each of these 
alternatives fails in some manner. 

The portion of the neutron star surface that is exposed to the hot fireball will drive 
a super-Eddington wind as the fireball begins to contract (TD95). This could force a 
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steady increase in the density of ions and electrons, supported above the surface by the 
photon pressure. At high enough temperatures the radiative flux from the surface of 
the fireball is fixed self-consistently by the opacity of this suspended matter (TD95): 
\asBT^ff - FeAd[<j{Teff,B)/aT\~''eion{i/Rm)~^- Here, Fedd is the Eddington flux 
calculated with the Thomson cross-section; a{Teff,B)/aT = 5'7i'^[kT/{heB/mec)]'^ is the 
scattering opacity in the strong background magnetic field; Eion is the ion density compared 
with the maximum that can be supported against gravity; and i is the thickness of the 
ion-electron layer at the photosphere. As the fireball contracts and i increases (due to the 
increasing exposure of heated crust) the radiative flux decreases as and the effective 
temperature decreases as Tg// oc Although this variant offers a nice explanation 

for the decrease in radiative flux and temperature during the smooth ~ 40 s tail, it has 
difficulty accomodating the contraction of the photosphere that is required by the transition 
to large-amplitude pulsations. 

Neutrino cooling does not change the lightcurve of a homogeneous fireball in the 
manner needed to explain the ~ 40-s smooth tail (Figs. ^b]) . The fireball could, of 



course, have two components, one of which has a much higher temperature Thigh and emits 
neutrinos much more rapidly than the cooler component. Then Thigh t~^^^ and the 
radiative flux from the contracting surface of the fireball decreases as(^ F oc Tf^lgh oc 
However, it is difficult to understand why the energy radiated from the higher temperature 
component should be much smaller (by a factor ~ 20), even while it covers a much larger 
volume (so as to explain the absence of large-amplitude pulsations). 



8. The Four-Peaked Repetitive Pattern: Collimated X-ray Jets 

A dramatic four-peaked pattern emerged in the August 27 light curve at ~ 40 s 
following the burst trigger. This pattern repeated with the 5.16-s rotation period of the 
neutron star. At times, the flux varied by more than an order of magnitude from peak to 
trough. However, there is no significant break in the pulse-averaged light curve during the 
appearence of the four-peaked pattern (Fig. 4 of Paper I), which indicates that the beamed 
flux is redistributed over the pulse period. 

The phase stability of the X-ray jets suggests that they are tied to surface features on 
the neutron star. Indeed, it has previously been argued that the radiative flux out of a 
trapped fireball could become significantly collimated along the magnetic field lines that 
open out to a few neutron star radii (TD95). Figure |^ outlines the basic geometry. A 



^^This applies when the higher Landau levels are populated; eq. (89) of TD95. 



-47- 



10^ — 10^^ G magnetic field is certainly strong enough: its pressure exceeds the radiative 
momentum flux by some 9-10 orders of magnitude. In this picture, the approximate ~ 1-s 
periodicity apparent in the four-peaked pattern is a chance byproduct of the location of the 
outburst, and the multipolar structure of the neutron star's surface magnetic field (Paper I; 
Fig. 1^). The difference in the number of sub-pulses observed within each rotation during 
the August 27 and March 5 giant flares, can then be ascribed to a difference in the number 
of xray jets. The most plausible geometry for each jet is a fan beam, which is swept past 
the line of sight once or twice by the rotation of the star. Figure ^ illustrates one possible 
geometry for each event, where each fan beam is observed twice as two separate sub-pulses. 

The large disparity between the scattering cross sections of the E-mode and 0-mode is 
the underlying reason for collimation. In addition, the scattering opacity of the E-mode rises 
rapidly with radius, as oc B^^ oc i?^ in a dipole geometry. As a result, radiative transport 
across the magnetic field lines is concentrated close to the neutron star surface (TD95). 
The trapped fireball heats a thin outer skin of the neutron star crust, and as the fireball 
contracts the cooling flux of X-rays drives matter off its surface (Ibrahim et al. 2000). This 
ablated material is easily suspended by Compton scattering]^ in the magnetosphere, where 
it remains confined. 

The rapid growth of the E-mode opacity then provides a mechanism for self-collimation: 
the E-radiation can escape only by pushing the suspended matter to the side (Fig. 
The ion-electron photosphere of the fireball is congruent with a set of magnetic field 
lines, as the result of pressure gradient forces along the field. (The cooling time exceeds 
the sound-crossing time of the fireball by a factor of a million.) This means that the 
collimation occurs primarily along magnetic field lines that open out beyond the electron-ion 
photosphere. 

The width of each X-ray jet depends on the amount of matter advected with the 
photons, and can be estimated as follows. A significant fraction of the E-mode flux near 
the E-mode photosphere is converted to the 0-mode by photon splitting (TD95). Mode 
changing also occurs via non-resonant Compton scattering near the E-mode photosphere 
( [Miller 1995| ; TD95), as well as by resonant Compton scattering near the ion cyclotron line 
(Thompson 2000). The 0-mode photons have scattering cross-sections near Thomson; i.e., 
the scattering of 0-mode photons is not significantly suppressed by the strong magnetic 
field. Thus, the energy flux injected into the 0-mode near the neutron star surface is 
tremendously super-Eddington, both due to the large luminosity Lq ~ \Lx and the small 



^^By non-resonant scatterineg off the 0-mode, or by resonant scattering at the ion cyclotron fundamental 
(which lies in the X-ray range in ~ 10^'' — 10^^ G magnetic fields; Ibrahim et al. 2001). 
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beaming angle Af2jet: 
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10^2 ergs s-i / V 47r 
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The 0-mode flows hydrodynamically along the magnetic field even in the presence of a tiny 
amount of matter, which can generate a large scattering depth along the magnetic field: 
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The resulting radiatively-driven outfiow is probably relativistic (§6.4 of TD95). Indeed, 
it quickly becomes relativistic if the photons dominate the specific heat and their expansion 
is adiabiatic. The conserved enthalpy per photon can be expressed in terms of the (bulk 
frame) photon density and the bulk Lorentz factor 7 as 7(^0)71^^ ~ kT{R^^)] and the 
rate at which photons are carried through the jet is 



= (7n^) A^et R^c 



Combining these two expressions gives 
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A^et(i?Ns)^: 
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which is 7 ~ -R/i?NS in spherical geometry but 
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in an outfiow that is channeled along a dipolar magnetic field]^. 



Photons travelling along field lines, with the angle between wavevector k and B 
satisfying 6'kB < {hu / meC^Y/'^{B / Bq-ed)'^/"^ ., experience strong suppression of photon 
scattering even in the O mode (e.g. Meszaros 1992). Inside this wavevector cone, the O and 
E photon eigenstates are nearly circularly polarized, with opposing helicities, rather than 
linearly polarized as at more oblique propagation angles. However, this cone is very narrow 
near the surface of the star: ^^kB < 0.05i?]~5^^^ (/;.a;/30keV)^''^ radians. In the zone where 
the jet outfiow is accelerated, just above the stellar surface, the 0-mode specific intensity is 
nearly isotropic, and almost all photons propagate outside the cone. 



13 



In this geometry, the cross-sectional area of a flux bundle increases as ^ 
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As the matter accelerates and expands, the optical depth through it drops off. Let us 
estimate the radius Rr where the matter and photons decouple, under the assumption that 
^kB ^ 1 at this radius in the bulk frame. This decoupling radius must be smaller than the 
maximum radius of the confining magnetic field lines {9^2 -^ns in a dipole geometry; Fig. 

< = (98) 



To reach a radial path of polar angle and escape, the advected photons must diffuse 
through an angle ~ \6 from the direction of the bulk streaming (which is tangent to the 
dipolar magnetic field). Thus the photons must diffuse a transverse distance Ai?^ ~ \9R, 
which corresponds to an optical depth 

Ti_ ~ ^^^^t^arAi?^. (99) 

7 

The diffusion time t^R^/c must be shorter than the radial flow time Rj^c, and we deduce 

n^{R) i^-e] otR < 1. (100) 



(101) 



,2 

Making use of eqs. (p4|) and (|99|) , this becomes 

Rr Ai7jet(-RNs)''"||(-RNS 



Combining this result with eq. (p8|), we find that the narrowness of the jet at its base 
increases increases with the optical depth: 

A^et(i?Ns) ^ 1 (^Q2) 



47r 2Ty\R 



NSj 



Further coUimation can occur when the escaping X-rays cross the surface of the electron 
cyclotron resonance, which sits at a distance R/Rns = 6.7(/itu/40 keV)^^/^ (Sns/IO^^ G)^/^ 
from a neutron star with polar field -Bns- During an SGR outburst, the matter suspended 
and confined at this radius can easily generate a large Thomson optical depth. This 
material is constrained to move along the magnetic field (the magnetic pressure greatly 
exceeds the radiative energy flux even at such a distance), and so its angular distribution 
will reflect the multipolar structure of the magnetic field near the stellar surface. However, 
in the persistent emission, which has a flux well below Eddington, the X-ray pulse profile 
will be strongly modified by resonant cyclotron scattering off magnetospheric currents (at 
either the electron or ion resonances: Thompson et al. 2001). 
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9. Conclusions 

The August 27 giant flare provides a Rosetta stone for SGR 1900+14, just as the 
March 5 giant outburst did for SGR 0526-66. These two remarkable events share very 
similar peak luminosities, energies, and morphologies, which suggests in turn that these two 
SGR sources are fundamentally alike. 

In the context of the magnetar model, we have discussed the mechanism by which such 
a giant flare may be triggered, and the physical processes operating in each of the three 
principal phases of the August 27 flare: the initial ~ 0.4 sec hard spike; the intermediate 
smooth ~ 40 sec tail, and the flnal phase of large-amplitude pulsations that displays a 
striking four-peaked pattern. In the process, we have reflned physical arguments that point 
to an external (dipole) magnetic fleld stronger than ~ 10^^ G in SGR 1900+14, and an even 
stronger ~ 10^^ G internal (toroidal) fleld which is the basic energy source for repeated flare 
activity. 

9.1. Physics of giant flares 

1. In the flare mechanism proposed here, most of the potential energy that powers 
a giant flare is stored before the event in the magnetic fleld of the deep crust and liquid 
stellar interior (TD95). This slowly-evolving fleld strains the crust from below. The brittle 
crust acts as a gate for the catastrophic release of energy. The elastic energy released in the 
crust is 5ma//er than the available magnetic energy by a factor ~ B'^/Awfi, where B is the 
fleld in the stellar interior and /i the crustal shear modulus, and by appropriate geometrical 
factors (§13.11). 

2. When the internal magnetic fleld is strongly wound, its stored energy can be rapidly 
communicated to the stellar exterior via a propagating fracture, involving a rotational 
deformation of a patch of the crust (§ p.2|) . Such a motion of the neutron star crust can 
create tangential discontinuities in the magnetic fleld, and thereby induce dissipation in 
three distinct zones. First, it induces strong magnetic shear in parts of the magnetosphere, 
which can rapidly damp through reconnection and conversion to high frequency Alfven 
waves; second, a torsional oscillation of the magnetized core is excited (along with an 
accompanying toroidal deformations of the crust which involve less energy); and third, 
static current sheets are excited deep in the crust. 

3. The damping of an internal shear mode by a flux of Alfven waves into the 
magnetosphere has been further quantifled. In the case of a large-scale mode with harmonic 
^ ^ R^s/'^Rfi (eq. |15), the equilibrium energy stored in a trapped magnetospheric Alfven 
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mode is greatly suppressed with respect to the exciting internal mode (eq. [0). This 
implies strong lower bounds on the poloidal magnetic field {B ^ 10^^ G) and internal mode 
energy 5-Eeias ^ 10'^^ ergs) needed to power the extreme peak luminosity of the initial hard 
spike ( §5.2| ). We also make a critical comparison of our results with the vibrating neutron 
star model of Ramaty et al. (1980), and review why an internal /- or p-mode couples even 
more weakly to the magnetosphere than an internal shear mode. 

4. The contracting photosphere of a very hot (T ~ 1 MeV) confined fireball provides 
an excellent fit to the envelope of the 1998 August 27 giant fiare after ~ 40s, and accounts 
for the rapid final drop in fiux (Paper I). We show in §^ that the fireball probably formed 
hotter in its center than at its edge, because the observed fiux diminished faster than it 
would for a homogeneous, spherical fireball. 



5. The light curve appears not to be perturbed by neutrino cooling (Figs. |4^, ^), 
which implies an upper bound to the fireball temperature of T ^ Tmax = 0.8 MeV. This 
in turn provides a conservative lower bound to the volume and magnetic moment of the 
confining magnetic field. The August 27 outburst alone implies a strict lower bound 

fi > 2.4 X 10^"*^ ergs) {T/T^g^x}^'^ G-cm'^. However, two considerations suggest that 

the net dipole moment is several times larger: the inferred number ~ 10^ of giant outbursts 
over the history of the source; and the absence of a measureable change at the ~ 10 percent 
level in the long-term spindown of SGR 1900+14 (Woods et al. 1999; Thompson et al. 
2000). 

6. The smooth ~ 40-sec tail of the 27 August fiare is somewhat harder than the 
ensuing large-amplitude pulsations. During this intermediate phase, the modest fiux 
variations do not repeat coherently with rotational phase; and the X-ray fiux is significantly 
higher than predicted by the trapped fireball model which fits observations after ~ 40 sec. 
In Paper I, we ascribed this excess hard fiux to Compton heating by an extended pair 
photosphere, driven by a persistent seismic excitation of the neutron star. In this paper, 
we have quantified the behavior of a pair corona in a magnetic field stronger than -Bqed- 
Cooling occurs primarily by Compton heating of the 0-mode. We defined a critical coronal 
luminosity (eqs. |82| and ^) below which a steady balance between electrostatic heating 
and diffusive radiative cooling is possible. This luminosity is 0(10^^ ergs s~^), which is 
comparable to the observed output from SGR 1900+14 during the smooth tail. A magnetar 
fiare could inject energy into the magnetosphere at a much higher rate than this critical 
level. Indeed, the fiux of the prompt, hard spike exceeds this minimum luminosity, giving 
direct evidence that energy was liberated fast enough to form a trapped fireball at the onset 
of the August 27 fiare. 



7. The best fit black body temperature is stable during the period of large-amplitude 
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pulsations (Mazets et al. 1999; Paper I) and agrees well with the value fcTgp ~ 11 keV at 
which photon splitting freezes out (eq. [180] in TD95). 

8. We propose that the large- amplitude pulsations are due to a collimated flux of 
X-rays from the base of the trapped fireball, moving along extended magnetic field lines 
(see also TD95; Feroci et al. 2001). CoUimation of the E-mode is provided by the rapid 
increase in its scattering opacity with distance from the stellar surface, (7e{B) oc B^^ oc . 
The 0-mode flows hydrodynamically even in the presence of a small flux of advected ions 
(and neutralizing electrons), Mc^/Lq <^ (Lo/Ledd)~H^^Ns/c^)~^- Energy can be released 
in both polarization modes at comparable rates. The width of the X-ray 'jet' can be related 
directly to the flux of advected matter. The collimation becomes finer as the matter flux 
increases, because the radiation-hydrodynamical flow must extend to larger radius before 
the photons and matter can decouple. 

9. We have quantified some physical mechanisms that could generate the persistent, 
hard spectral component of the August 27 flare. In particular, the burst spectrum is harder 
during dips, which points to acceleration of non-thermal particles or direct Comptonization 
by large-amplitude Alfven waves in an extended corona. 



9.2. Evidence for Magnetars: Field Strength Estimates 

One of the principal goals of this paper is to set more model-independent (lower) 
bounds on the magnetic field in the flaring SGR sources. Several observational properties of 
these sources directly require magnetic fields stronger than ~ 10^^ G (the dipole component) 
to ~ 10^^ G (the internal toroidal field, and higher multipoles). 

1. The estimated output of ~ 10^"^ ergs in giant flares over the active history of a 
flaring source corresponds to a r.m.s. magnetic field stronger than ~ 10^^ G (§p.2|). 



2. The observed lightcurve of the August 27 flare can be well fit by the cooling of 
a 'trapped fireball' (Feroci et al. 2001). The confinement of a substantial fraction of the 
outburst energy [~ (1 — 3) x 10"^^ ergs] in a pair-photon plasma close to the source, which 
cools gradually over ~ 300 seconds, implies a lower bound ~ 10^^ G to the external dipole 



field (TD95). This argument is refined in § |5.1| , taking into account the fit of the August 



27 flare light curve to the contracting surface of a trapped fireball (Paper I), and the 
restrictions on the size of the fireball from bulk neutrino-pair cooling. In addition, the 
intermediate ~ 40 smooth tail in the August 27 flare provides evidence that the seismic 
output of the source can exceed ~ 10'^ Ledd for a limited time. 
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3. The extreme peak luminosity of the initial spike (3 — 10 x lO^Ledd^ Fenimore et al. 
1996; Mazets et al. 1999) can be powered through a sudden readjustment of a magnetic 
field stronger than ~ 10^"^ G (§p.2D, especially if this readjustment involved a fracture of 
the neutron star crust and the formation of regions of strong magnetic shear. Such a violent 
event would also deposit energy into internal shear modes (including a torsional Alfven 
wave in the liquid core and standing shear waves in the rigid crust). We have considered 
separately the damping of these internal oscillations via a coupling to the magnetophere. 
Even under the most optimistic assumptions about the rate at which an external torsional 
mode is converted to radiation, the observed peak luminosity of the giant flares can barely 
be supplied if the surface field is ~ 10^^ G, and cannot be supplied if it is much weaker. 
The luminosity of the pulsating tail, L/Lgdd ^ 10^, is consistent with the suppression of 
the electron scattering opacity of the E-mode radiation near the strongly magnetic surface 
of a trapped fireball. In the context of the SGR sources, this effect was first discussed by 
Paczyhski (1992), under the assumption that the radiation is released from the cooling 
surface of a magnetar. However, the physical effect is in fact much cleaner if radiative 
transport occurs across the confining magnetic field lines of a trapped fireball. 

It is also worth re-evaluating other arguments for strong magnetic field, in the light of 
recent observations: 

4. The long 8-s spin period of SGR 0526-66, combined with its association with 
the young LMC supernova remnant N49 (e.g. Cline 1982), suggests that this source is 
rapidly spinning down. Spindown from a much shorter period, in the age to526-66 of N49, 
corresponds to a polar dipole field of ~ 6 x 10^"^ (%26-66/10^ yr)~^^^ G if the classical 
magnetic dipole formula applies (DT92). Although rotational modulation of the persistent 
emission of SGR 0526-66 has not yet been detected with compelling statistical confidence 
(Kulkarni et al. 2000), two other sources SGR 1806-20 and SGR 1900+14 are observed to 
spin down rapidly (Kouveliotou et al. 1998b, 1999). 

The observed P variability of these stars (Woods et al. 1999; 2000) may be surprising 
at first sight, if they are identified as isolated, non-accreting neutron stars. However, as 
we have argued in detail in this paper, the bright X-ray outbursts provide direct evidence 
for sudden deformations of the external magnetic field in the SGR sources. The rate of 
spindown of an isolated neutron star is controlled by the electrical current flowing across 
its speed-of-light cylinder. Variations in the external magnetic field will, therefore, cause a 
modulation of the spindown torque - either through a large-scale twisting of the external 
magnetic field (Thompson et al. 2001); or through a magnetically-powered wind (Thompson 
& Blaes 1988; Thompson et al. 2000). The wind model requires dipole fields ^ 1 x 10^^ G 
for the SGRs with measured spindown, unless the energy in a wind exceeds the observed 
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X-ray output by more than an order of magnitude. 

Note that several radiopulsars with spindown fields B^ipoie > -Bqed have recently 
been found. The present record value of B^ipoie inferred from P and P in a radiopulsar is 
1.1 X 10^^ G (polar field) for PSR 1814-1744 (Camilo et al. 2000). This means that the 
minimum magnetic moment which we have inferred from the flare physics is near the upper 

end of the range measured in radio pulsars. The dipole field of PSR 1814-1744 is a factor 
10 — 20 weaker than that inferred from the spindown of SGRs 1806-20 and 1900+14 if the 
magnetic dipole formula applies; but can plausibly be reduced by a factor of 3 — 10 in the 
presence of persistent currents (Kouvcliotou et al. 1998b, 1999; Harding, Contopoulos, & 
Kazanas 1999; Thompson et al. 2000, 2001). 

On the other hand, PSR 1814-1744 sits close in the P — P plane to the magnetar 
candidate AXP 1E2259+586. Nonthclcss, the spindown age of IE 2259+586 exceeds by 
a factor ~ 10 — 30 the age of the surrounding supernova remnant CTB 109 (e.g. Kaspi, 
Chakrabarty, & Steinberger 1999). Given the rapid spindown observed in the other SGR 
and AXP sources, this suggests that the spindown torque of IE 2259+586 has decayed by 
a factor ~ 10 from its historic average. Note that the interior, toroidal magnetic fields of 
radiopulsars could be much weaker than those of SGRs and AXPs (TD93), which would 
provide an explanation for why PSR 1814-1744 is a much weaker X-ray source than the 
AXPs (Pivovaroff, Kaspi, & Gotthelf 2000). 

5. The absence of a significant perturbation to the long-term spindown rate of SGR 
1900+14 in the few months following the August 27 event, imphes only a small change 
in the external dipole field coinciding with the release of ~ 10^^ ergs (Woods et al. 1999; 
Thompson et al. 2000). 

6. The detection of extended afterglow from the heated surface of a Soft Gamma 
Repeater is expected following exposure to a trapped fireball. The amount of fireball energy 
absorbed by the crust increases linearly with the surface B-field at large magnetic fiux 
densities (TD95). This effect has been observed in some short SGR bursts (Strohmayer & 
Ibrahim 1998), and has been related to an extended, faint oscillatory tail that was observed 
in an August 29 outburst from SGR 1900+14 (two days after the August 27 giant fiare; 
Ibrahim et al. 2000). 
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Table 1. Fireball Indices 
Dimensions Energy Density Flux Fireball Structure 
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Note. — The last two columns give the values of fireball 
structure indices needed to explain the observed time-dependence 
of the August 27 event = 3). The symbol * means that the 
observed time-evolution is insensitive to this index, i.e., any value 
would work. 
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Fig. 1. — An almost spherical star is threaded by a uniform poloidal magnetic field. This 
field threads both the hquid core and the solid, outer crust of the star. In the core, the field 
is twisted, creating a toroidal component. 
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Fig. 2a. — A twisted bundle of magnetic flux threads a conducting medium, which is 
stratified along planes which run perpendicular to the direction of gravity. The flux bundle 
is tilted by angle a with respect to gravity; it may be anchored from above and below. 
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Fig. 2b. — When the flux bundle runs perpendicular to the direction of stratification, a 
localized twist is prevented from spreading out along the field through purely hydromagnetic 
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FIREBALL BOUNDARY 




a SURFACE b 

Fig. 3. — A trapped fireball releases energy as its cool outer surface contracts. We illustrate 
two possible geometries (approximately spherical and cylindrical) for the fireball. The 
measured fluence of the extended tail of a giant flare, combined with the absence of a 
measureable perturbation to the light curve from pair neutrino cooling, sets a lower bound 
to the dipole moment of the confining magnetic field. 
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Fig. 4a. — Surface X-ray luminosity of a homogeneous, trapped fireball, whose radiative 
area scales with the fireball volume as ~ ('^ = 3. lower curves) and A ~ V^-^ (x = 1; 
upper curves). Within each set of curves, the solid curve includes no neutrino cooling. The 
long-dashed and short-dashed curves correspond to uniform pair neutrino cooling, with an 
initial cooling rate 1 and 10 times the initial surface X-ray luminosity. 
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Fig. 4b. — Surface X-ray luminosity of a spherical trapped fireball {Dc — 2), showing the 
effects of neutrino cooling on the shape of the hghtcurve. The bold curve (x = 3) closely 
approximates the observed lightcurve of the August 27 flare. In a homogeneous magnetic 
field, this fireball index corresponds to a mild temperature gradient 1 = —\- The short- 
dashed curve describes the effect of pair-neutrino cooling on the X-ray lightcurve, with the 
neutrino luminosity initially 100 times the cooling X-ray luminosity. The same relative 
neutrino cooling rate, but a different temperature structure (7 = 0), brings closer agreement 
with the observed light curve, but there remains some difference in shape (long-dashed curve). 



-68- 



Fig. 5. — A persistent seismic excitation of the neutron star crust will in turn excite sheared 
(current-carrying) Alfven waves on extended magnetic field lines in the star's magnetosphere. 
Non-linear couplings between these trapped Alfven modes generate a turbulent cascade that 
creates a hot corona. At the high luminosity (and inferred compactness) of the August 
27 smooth tail, this corona consists of electron-positron pairs with only a slightly higher 
temperature than the X-ray photons (which are assumed to have a Wien distribution). The 
dominant source of fresh Compton seeds for this Corona could either be provided by an 
interior trapped fireball (as depicted); or, alternatively, by photon splitting. 
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Fig. 6. — Equilibrium temperature of a electron-positron corona (of size R), under two 
different assumptions about the mechanisms of photon creation and radiative cooling. The 
top and bottom curves correspond to an external source of (E-mode) photons (eq. ||7^). The 
second set of curves corresponds to an internal photon source, which is balanced by outward 
diffusion through the 0-mode (eq. ^9|). The maximum coronal luminosity, above which no 
stable balance between heating and cooling is possible, lies close to ~ 10"^^ erg s~^ (eqs. |^ 
and [|9|). 
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Fig. 7. — A trapped fireball releases energy as its cool outer surface contracts. The radiative 
fiux across the fireball surface is concentrated close to the surface of the neutron star, where 
the E^modc scattering opacity is most strongly suppressed by the intense magnetic field. 
A collimated flux of both the 0-mode and E-mode radiation can be generated along more 
extended magnetic fleld lines. The 0-mode couples tightly to the electrons near the stellar 
surface, even if a small fraction of the hvDer-Eddineton radiative flux is carried bv matter: 
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Fig. 7. — A geometry is suggested both for the August 27 giant flare from SGR 1900+14 
(left; four sub-pulses) as well as the original March 5 flare from SGR 0525-66 (right; two 
sub- pulses). Each large arrow denotes a fan beam of X-rays and outflowing, relativistic 
matter. The angular velocity of each source is assumed to be oriented so that each fan beam 
is observed twice per rotation, yielding two sub-pulses per beam. In this model, the burst 
light curve is sensitive to the configuration of the particular subset of the closed magnetic 
field lines which confine the hot e^-photon plasma. Thus, the presence of higher multipoles 
in SGR 0526-66 is not excluded by the observation of two sub-pulses within the tail of the 
March 5 flare. 



